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(1)  Development  of  a  new  propagation  prediction  model 
applicable  at  frequencies  up  to  5  GHz. 

(2)  Development  of  a  modem  prediction  model  for  cal** 
culating  the  bit  error  rate  and  fade  outage  on  a 
fading  multipath  channel  for  the  MD-918/GRC  modem. 

(3)  A  trade-off  analysis  to  determine  the  optimum 
angle  diversity  system. 

(4)  Design  and  development  of  a  4. 4-5.0  GHz  (C-Band) 
angle  diversity  feedhorn  which  gives  a  beam 
separation  within  about  1/2  dB  of  optimum. 

(5)  Theoretical  and  empirical  evidence  that  angle 
diversity  has  a  performance  advantage  over 
frequency  diversity. 

(6)  Discovery  of  reduced  correlation  effects  in  angle 
diversity  when  wideband  signals  are  employed. 

(7)  Development  of  a  methodology  and  a  statistical 
base  for  development  of  future  angle  diversity 
systems .  |U 

This  effort  impacts  tne  development  of  future  digital  tropo- 
scatter  circuits  by  establishing  a  new  technology  which  will 
allow  conversion  from  frequency  to  angle  diversity  with  a 
resulting  halving  of  frequency  requirements  and  a  positive 
system  gain  in  almost  ail  applications.  The  angle  and  fre¬ 
quency  diversity  systems  have  virtually  the  same  level  of 
complexity. 
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SECTION  1 


INTRODUCTION 

The  Defense  Communications  Agency  (DC A)  and  the  military 
departments  are  upgrading  the  troposcatter  links  in  the  Defense 
Communications  System  (DCS) .  Present  efforts  are  focused  on 
DCS  links  in  Europe.  NATO  is  also  contemplating  an  upgrade  of 
the  troposcatter  links  in  its  ACE-High  system.  Current  plan¬ 
ning  for  DCS  troposcatter  links  calls  for  all  digital  traffic 
to  meet  the  criteria  set  forth  in  the  Defense  Communications 
Engineering  Center's  Technical  Report  No.  12-76,  "DCS  Digital 
Transmission  System  Performance".  The  upgrade  effort  is  com¬ 
plicated  by  a  need  to  reduce  frequency  assignments  and  con¬ 
serve  spectrum  bandwidth. 

In  response  to  DCA  tasking  for  a  solution  to  the  above 
problem  areas,  the  Department  of  Army  has  sponsored  an  in¬ 
vestigation  of  angle  diversity  as  a  promising  technique  to 
combat  short  term  fading,  and  alleviate  aperture-to-medium 
coupling  loss  without  expanding  and  possibly  reducing  current 
usage  of  the  frequency  spectrum.  The  effort  consisted  of 
three  phases:  the  study  phase,  the  equipment  development 
phase  and  the  link  test.  In  brief,  these  phases  consisted 
of  the  following: 

•  During  the  study  phase,  comprehensive  models  of 
the  troposcatter  propagation  and  digital  modem 
characteristics  were  modeled.  The  propagation 
model  emphasized  the  turbulent  scatter  mechanism 
as  the  limiting  mode,  particularly  for  S  and  C 
band  (2. 0-5.0  GHz)  transmissions.  This  model  was 
written  in  terms  of  spectrum  slope,  refractive  in¬ 
dex  variance,  and  scale  of  turbulence.  In  the  past, 
aperture-to-medium  coupling  loss  has  been  included 


in  total  path  loss  calculations  as  a  separate  parameter, 

i.e.  the  total  path  loss  is  the  basic  path  loss 

with  omnidirectional  antennas,  the  antenna  gains, 

and  the  coupling  loss.  The  usefulness  of  this 

concept  is  questionable,  since  it  uses  properties 

of  the  atmospheric  structure  outside  the  actual 

antenna  beams.  The  approach  developed  in  this 

study  calculated  the  total  path  loss  directly  by 

integrating  over  the  common  volume.  The  digital 

modem  model  was  developed  to  predict  bit  error 

rate  (BER)  statistics  for  the  MD-918  troposcatter 

modem.  This  modem  employs  a  decision-feedback 

equalizer  as  the  means  of  combining  diversity 

channels  and  minimizing  intersymbol  interference 

effects.  The  MD-918  modem  has  been  tested  in  DCS 

troposcatter  configurations  in  both 

the  United  States  and  in  Europe.  The  propagation 

and  digital  modem  models  were  used  to  predict 

the  performance  of  angle  diversity  in  replacing 

space  or  frequency  diversity  and  in  augmenting 

either  or  both  for  systems  providing  up  to  eighth 

order  diversity. 

Equipment  Development  Phase:  The  angle  diversity 
design  was  selected  as  a  vertical  splay  of  two 
beams  spaced  by  approximately  1  beamwidth.  A 
C-Band  angle  diversity  feed  horn  was  developed 
with  capability  of  duplex  operation  with  up  to 
10  Kw  transmitters  and  with  dual  polarization  ports 
in  both  horn  branches.  A  digital  traffic  capa¬ 
bility  was  realized  through  the  use  of  MD-918 
digital  tropo  modems.  A  predetection  combiner 
(PDC)  was  developed  to  implement  up  to  8th  order 
diversity  with  the  MD-918  for  combining  the  signals 
of  the  elevated  (i.e.,  the  angle  diversity)  beams. 


•  Link  Test  Phase:  The  test  program  consisted  of 

an  extensive  series  of  propagation  and  modem  performance 
tests  from  October  1977  to  May  1978  followed  by 
an  effort  to  collect  long  term  propagation 
data.  This  second  effort  continued  until  October 
1979.  The  test  bed  for  these  experiments  was  the 
RADC  168  mile  C-band  troposcatter  system  because 
Youngstown  and  Verona,  N.Y.  Specific  parameters 
measured  were  hourly  path  loss  distributions, 
long  term  distributions  of  hourly  median  path 
loss,  multipath  profiles  and  rms  multipath  spread, 
short  term  correlation  between  angle  diversity 
paths,  long  term  decorrelation  of  angle  diversity 
hourly  medians,  angle  diversity  squint  loss  sta¬ 
tistics,  and  BER  statistics  at  6.3  and  12.6  Mb/s. 

The  major  results  of  this  effort  are  best  summarized  with 
respect  to  the  program  phase  in  which  it  occurred. 

a.  Study  Phase  Results 

1.  A  new  troposcatter  propagation  prediction  model 
applicable  at  frequencies  up  to  5  GHz  was 
developed.  The  present  NBS  Tech  Note  101  method 
for  propagation  prediction  is  a  semi-empirical 
method  limited  to  carrier  frequencies  of  1  GHz 
or  less. 

2.  A  method  for  predicting  BER  outage  rate  was 
developed  for  the  MD-918  digital  troposcatter 
modem.  Previous  performance  calculations  were 
based  on  an  average  probability  of  error  con¬ 
cept  which  is  a  poor  measure  for  digitized 
speech.  DCS  standards  are  now  specified  in 
terms  of  Bit  Error  outage  rate. 
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3.  A  trade-off  analysis  was  completed  to  determine 

the  best  angle  diversity  configuration.  The  op¬ 
timized  system  consists  of  a  main  beam  transmit-  A 

ter  for  each  antenna  and  two  vertically  splayed 

beams  for  each  antenna  at  the  receiver.  The  an¬ 
tennas  are  aimed  with  boresight  slightly  above 
the  horizon  as  with  conventional  dual  space/dual 
frequency  diversity  (2S/2F)  systems.  A  vertical 
splay  of  about  3/4  beamwidth  is  desirable  al¬ 
though  physical  limitations  in  feedhorn  design 
usually  result  in  a  somewhat  larger  splay.  The 
dual  space/dual  angle  (2S/2A)  diversity  system  interfaces 
with  the  digital  modem  in  the  same  manner  as  the 
conventional  2S/2F  system,  i.e.,  no  modem  changes 
are  required. 

4.  A  theoretical  comparison  of  2S/2F  and  2S/2A  di¬ 
versity  configurations  showed  the  angle  systems 
to  be  generally  better  than  the  frequency  diver¬ 
sity  systems.  The  importance  of  this  result  is 
reflected  in  the  savings  of  two  bandwidth  alloca¬ 
tions  per  link  in  a  conversion  from  frequency  to 
angle  diversity. 

b.  Equipment  Development  Phase  Results 

1.  A  C-band  (4.4  -  5.0  GHz)  dual  polarized  feedhorn 
with  a  beam  separation  of  1.3  beamwidths  was  de¬ 
veloped.  This  beam  separation  results  in  a  sys¬ 
tem  loss  of  about  1/2  dB  relative  to  the  optimum 
separation  of  3/4  beamwidth.  The  insertion  loss 
due  to  replacing  the  original  horns  with  the  an¬ 
gle  diversity  horns  was  on  the  order  of  0.1  dB. 

2.  An  extension  of  the  MO-918  digital  modem  was  de¬ 
veloped  to  accommodate  up  to  eight  diversity  in¬ 
puts.  The  eighth  order  diversity  system  results 

when  2S/2F  is  augmented  with  angle  diversity  to  * 

produce  2S/2F/2A. 
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g.  Link  Test  Results 


1.  Angle  diversity  was  shown  to  have  better  performance 
than  frequency  diversity. 

2.  Eighth  order  diversity  offers  only  a  small 
improvement  over  fourth  order  systems. 

3.  The  short  term  correlation  between  angle  diver¬ 
sity  beams  has  an  insignificant  effect  on  per¬ 
formance  for  wide  bandwidth  (~5  MHz)  digital 
systems . 

4.  Elevated  beam  multipath  spread  exceeds  main 
beam  multipath  spread  by  about  a  third  or  less. 

The  distribution  slopes  of  the  multipath  spread 
are  about  the  same. 

5.  The  elevated  beam  frequently  has  a  stronger 
received  signal  than  the  main  beam.  Since  the 
main  beam  has  been  aimed  for  maximum  received 
power,  this  result  requires  an  inhomogeneous 
atmospheric  structure  with  respect  to  the  two 
common  volumes. 


d.  Conclusions 

These  results  impact  on  the  development  of  future  digital 
troposcatter  systems  in  a  number  of  ways.  Because  analog 
FM/FDM  troposcatter  systems  generally  require  less  bandwidth  than 
digital,  by  almost  a  factor  of  two,  the  conversion  to  digital  of  2S/2F 
X  systems  requires  twice  as  many  bandwidth  allocations.  If, 

however,  the  digital  conversion  is  accompanied  by  a  conversion 
of  the  frequency  diversity  to  angle  diversity,  the  same  diver¬ 
sity  performance  is  realized  but  with  the  same  total  bandwidth  re- 
•  quirements  as  the  analog  system  rather  than  twice  these 

requirements.  In  geographic  areas  such  as  Europe,  where  fre¬ 
quency  assignments  are  hard  to  obtain,  digital  conversion 
may  not  be  possible  without  use  of  angle  diversity.  The 
t  performance  advantage  of  angle  over  frequency  diversity  has 


* 
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not  been  recognized  in  the  past  because  of  the  emphasis  on  the 
squint  loss  and  correlation  of  the  elevated  beam.  Our  results 

*  ► 

show  that  in  wideband  digital  systems,  correlation  is  not  impor¬ 
tant  and  that  the  squint  loss  is  compensated  by  both  the  second 
transmitter  power  amplifier  which  becomes  available  when  a  fre¬ 
quency  to  angle  diversity  conversion  is  made  and  an  availability  gain 
of  angle  diversity  due  to  decorrelation  in  the  two  beams  common  vol¬ 
umes.  Two  power  amplifiers  are  used  in  DCS  communication  applications - 
to  insure  link  reliability  by  providing  redundancy. 

The  digital  conversion  of 

will  involve  some  links  with  carrier  frequencies  in  the  4.4 
to  5  GHz  band.  The  new  prediction  model  developed  under  this 
program  will  be  essential  for  determining  upgrade  recommendations 
and  calculating  network  performance.  The  use  of  previous  methods 
which  were  not  derived  for  this  frequency  range  may  have  led  to 
unsatisfactory  link  service. 

The  small  improvement  resulting  from  extension  to  eighth 
order  diversity  is  due  to  the  necessity  to  improve  the  outage 
rate  at  signal-to-noise  values  where  the  point  of  diminishing 
returns  with  respect  to  diversity  has  been  reached.  Thus 
systems  which  do  not  meet  outage  rate  standards  generally  ( 

require  brute  force  dB  improvement  such  as  larger  power 
amplifiers  or  antennas. 

The  absence  of  a  correlation  degradation  in  wideband 
digital  angle  diversity  systems  is  due  to  the  decorrelation  0| 

effects  of  multipath  on  the  two  beams.  Thus  a  correlation 
coefficient  of  0.6  measured  with  CW  signals  which  would  degrade  i 

a  narrowband  digital  system  by  1  dB  has  no  significant  effect 
on  a  wideband  system.  Thus  in  future  design  of  angle  diversity  4 

feedhoms  every  attempt  should  be  made  to  reduce  the  beam  j 

separation  and  hence  the  squint  loss.  Reduction  of  the  beam 
separation  should  not  be  at  the  expense  of  unacceptable  cou¬ 
pling  between  horns,  however.  1 
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The  multipath  spread  results  are  encouraging  m  that  on 
most  systems  the  present  MO-918  modem  can  accommodate  the  additional 
multipath  spread  associated  with  the  elevated  beam. 


e.  Report  Contents 

This  report  presents  the  significant  results  of  this 

program  in  the  following  manner.  Section  2  introduces  key 

troposcatter  concepts  and  develops  the  angle  diversity  design. 

The  propagation  model  is  presented  in  Section  3  along  with  a 

summary  of  the  computer  program  for  predicting  path  loss. 

Section  4  develops  the  performance  model  for  the  MD-918  digital 

troposcatter  modem.  In  Section  5  the  DCS  performance  standards 

are  reviewed  and  an  example  computation  of  link  performance  is 

made  and  compared  to  these  standards.  During  the  course  of 

this  program,  this  link  analysis  was  completed  for  all  identified 

DCS  and  NATO  links.  The  results  of  this  analysis  are  available 

in  a  previous  report The  last  section  addresses  the  question 

of  long  term  availability  and  the  potential  for  gain  in  an  angle 

diversity  system.  This  section  also  presents  key  test  results 

fl  21 

which  are  contained  in  greater  detail  in  1  .  Additional  details 

on  specific  phases  of  the  program  are* given  in  the  two  interim 
reports  *1,3'  which  are  available  from  the  Defense  Documenta¬ 

tion  Center. 
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SECTION  2 


ANGLE  DIVERSITY  IN  DIGITAL  TROPOSCATTER  SYSTEMS 

Studies  of  adaptive  antenna  control  techniques  under  this 
program  have  shown  that  a  small  gain  can  be  achieved  on  diffrac¬ 
tion  links  with  electro-mechanical  antenna  steering  and  that  a 
significant  savings  in  bandwidth  allocations  can  be  realized  on 
scatter  circuits  by  conversion  of  frequency  diversity  systems 
to  angle  diversity.  Because  of  the  importance  of  these  angle 
diversity  results,  this  report  emphasizes  the  troposcatter 
problem. 

Additional  angle-diversity  (AD)  channels  can  be  realized 

through  the  placement  of  extra  feedhorns  near  the  focal  point 

of  the  radio  antenna.  In  one  of  the  early  angle  diversity 

\2  11 

experiments,  Crawford  et  al. 1  *  J  confirmed  that  vertically 

displaced  feeds  provide  essentially  uncorrelated  channels  at 

both  460  and  4110  MHz.  Their  experiments  with  horizontally 

displaced  feeds  revealed  considerable  correlation  at  the  lower 

frequency,  but  again  essentially  uncorrelated  channels  at  the 

higher  frequency.  In  another  study  of  azimuthal  angle  diver- 

\2  21 

sity  Chisholm  et  al.  *  J  determined  correlation  coefficients 
at  two  separations  and  the  effective  diversity  gain.  Vogelman 
et  al.^2*3^  also  established  the  low  correlation  between  AD 
signals. 

T2  41 

Surenian1  *  J  considered  the  use  of  angle  diversity  at 
super  high  frequencies  (SHF) .  The  use  of  SHF  for  troposcatter 
is  limited  by  the  increased  path  loss.  Angle  diversity  prom¬ 
ised  a  workable  alternative  to  the  use  of  additional  or  larger 
antennas  to  combat  the  additional  path  loss.  His  experiments 
at  8  GHz  established  that  effective  diversity  can  be  achieved 
with  only  a  single  transmitter  beam  and  multiple  receiver 


beans.  The  correlation  between  beams  was  adequately  small 

(<0.6)  and  for  a  center  cluster  of  beams  only  a  1  or  2  dB 

median  difference  in  signal  strength  was  observed.  The  use 

of  angle  diversity  on  a  long  troposcatter  path  was  reported 
[2  51 

by  Travis  *  who  showed  that  there  was  substantial  signal 
strength  in  the  elevated  beams  and  that  the  fading  was 
essentially  uncorrelated.  The  correlation  coefficients  for 
horizontally  displaced  feeds,  however,  were  greater  than  0.6, 
reducing  the  usefulness  of  azimuthal  AD  on  such  a  system.  An 
angle  diversity  test  described  by  Monsen1  *  1  provided  data 
on  the  long-term  decorrelation  of  the  hourly  medians  and  found 
close  agreement  between  a  dual  frequency  and  triple  angle 
diversity  system.  The  relatively  large  squint  angle  (~1.6 
beamwidths)  in  this  experiment  made  dual  angle  poorer  than 
dual  frequency  diversity  for  nominal  antenna  pointing  angles, 
f  2  71 

Troitskiy1  *  J  also  remarked  on  the  long-term  hourly 
median  decorrelation  advantage  of  angle  diversity  and  provided 
correlation  results  of  a  long-term  test  program.  Angle  diver¬ 
sity  tests  reported  on  by  Gough,  et  al. showed  good 

\2  91 

agreement  with  early  theoretical  work  by  Koono ,  Hirai ,  et  al 1  J . 

Previous  work  on  angle  diversity,  however,  did  not 
appropriately  assess  the  impact  of  the  following  considerations: 

1.  Frequency  diversity  requires  two  power  amplifiers, 
whereas  a  theoretical  angle  diversity  requires  only 
one.  From  a  reliability  redundancy  viewpoint  or  an 
equal  total  transmit  power  argument,  the  systems 
should  be  compared  on  the  basis  of  two  power  ampli¬ 
fiers  in  each  system.  Any  operational  angle  diver¬ 
sity  system  would  require  a  second  power  amplifier 
as  a  reliability  spare.  This  fact  increases  the 
received  power  in  the  main  beam  by  3  dB  in  an  angle 
diversity  system.  This  3  dB  increase  offsets  up 
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to  a  6  dB*  squint  loss  in  the  combined  diversity 
output. 


& 


2.  The  effect  of  correlation  between  diversity  signals 
is  reduced  when  signals  are  wideband  with  respect  to 
the  reciprocal  of  the  multipath  spread  and  adaptive 
signal  processing  techniques  are  used. 

In  addition  to  these  factors,  the  potential  for  angle 
diversity  increases  when  one  considers  that  major  disadvantages 
of  an  analog  angle  diversity  system  either  disappear  or  become 
advantages  when  angle  diversity  is  applied  to  a  digital  system 
with  adaptive  signal  processing.  These  disadvantages  in  an 
analog  system  are  the  variation  in  delay  between  angle  diver¬ 
sity  beams  and  increased  multipath  spread  in  the  elevated  beam. 

In  analog  systems  these  effects  led  to  intermodulation  dis¬ 
tortion  which  sharply  limited  message  capacity.  In  digital 
systems,  delay  variation  can  be  compensated  for  by  an  adaptive 
equalizer  and  additional  multipath  can,  in  many  cases,  lead  to 
better  performance  as  a  result  of  added  implicit  diversity. 

This  result  assumes  that  an  equalizer  type  signal  processor 
with  a  delay  span  large  enough  to  accommodate  the  delay  vari¬ 
ation  and  multipath  spread  is  used.  In  almost  all  applications 
considered  in  this  study,  a  three  tap  equalizer  with  one-half 
QPSK  symbol  spacing  was  found  to  be  adequate. 

In  this  section  we  develop  the  major  design  considerations 
for  angle  diversity.  The  preceding  discussion  has  established 
the  importance  of  adaptive  signal  processing  in  the  angle  diversity 
application.  Thus,  before  proceeding  with  the  design  consider¬ 
ations,  a  review  of  adaptive  systems  on  digital  troposcatter  is 
provided.  This  review  addresses  in  particular  the  impact  of  these 
systems  on  analog  to  digital  conversion  of  troposcatter  links  and 
the  role  of  angle  diversity. 


t  *A  squint  loss  of  S  dB  results  in  a  S/2  dB  combining  loss  because 

the  error  rate  is  a  function  of  the  geometric  mean  of  the  signal- 
to-noise  ratios. 
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a.  Adaptive  Signal  Processing  for  Digital  Troposcatter  Systems 

Troposcatter  radio  transmission  was  discovered  after  World  * 

War  II  when  it  was  noted  that  microwave  signals  from  beyond  the 
horizon  radars  were  much  stronger  than  predicted  from  diffraction 
calculations  over  the  earth's  surface.  A  commonly  held  theory 
of  this  phenomenon  is  explained  by  Tatarski .  The  theory 
is  that  random  fluctuations  in  the  dielectric  constant  in  the 
troposphere  divert  some  small  fraction  of  the  impinging  energy 
back  to  the  receiver.  The  name  tropospheric  scatter  or  tropo¬ 
scatter  derives  from  this  concept  of  random  redirection  of  the 
incident  wave  by  the  troposphere.  Significant  scatter  returns 
occur  from  a  "common  volume"  defined  by  the  antenna  beam  patterns. 
Scatter  returns  from  different  points  within  the  common  volume 
have  different  path  delays ,  and,  when  these  signal  scatter  points 
are  separated  by  more  than  the  decorrelation  distance  of  the 
dielectric  constant  fluctuations,  the  scatter  returns  are  not 
correlated.  The  impulse  response  characterizing  this  channel 
has  a  time-varying  continuous  multipath  structure.  A  statis¬ 
tical  measure  of  the  width  of  the  troposcatter  channel  impulse 
response  is  the  multipath  delay  spread. 

Troposcatter  systems  have  been  widely  used  in  military 
applications  for  beyond  line-of-sight  communications  up  to 
about  600  miles.  The  frequency  range  for  this  application 
extends  from  about  400  to  5000  MHz.  Reliable  communications 
require  redundant  transmission  paths  provided  through  the  use  3 

of  multiple  frequencies,  antennas  separated  in  space,  or  scatter 
at  two  different  beam  angles.  The  several  redundant  paths  are 
referred  to  as  diversity  paths.  The  number  of  paths  is  termed 
the  order  of  diversity.  The  multipath  delay  spread  limits  the  9 

channel  capacity  which  can  be  achieved  in  present  analog  sys¬ 
tems.  Only  transmission  bandwidths  less  than  the  reciprocal 
of  this  multipath  delay  can  be  achieved.  Signals  of  larger 
bandwidths  become  distorted  due  to  the  multipath  dispersion.  % 
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In  FM  systems  this  dispersion  causes  intermodulation  noise 
after  detection. 


With  the  introduction  of  satellite  communication  systems, 
which  do  not  suffer  from  extensive  multipath  fading,  the  future 
of  troposcatter  systems  appeared  to  be  limited.  Economic  and 
security  factors  have  altered  this  assessment,  particularly  for 
digital  transmission.  With  digital  signal  formats,  adaptive 
methods  can  be  devised  to  measure  the  multipath  structure  and 
exploit  it  as  an  extra  form  of  diversity  to  improve  performance. 
Unlike  the  capacity. of  analog  systems,  the  capacity  of  digital 
systems  is  not  restricted  by  the  reciprocal  of  the  multipath 
width.  Also  from  a  network  viewpoint  in  a  digital  system, 
fades  in  tandem  links  do  not  have  a  cumulative  effect  because 
the  signal  can  be  regenerated  at  each  node. 

Applicability  of  adaptive  signal  processing  techniques 
is  critically  dependent  on  whether  the  rate  of  fading  is  slower 
than  the  rate  of  signaling.  For  transmission  rates  of  interest, 
troposcatter  radio  links  can  be  considered  to  be  slow  fading 
multipath  channels. 

(1)  The. Learning  and  Diversity  Constraints 

For  digital  communication  over  beyond-the-horizon  radio 
links,  an  attempt  is  made  to  maintain  transmission  linearity, 
i.e.,  the  receiver  output  is  a  linear  superposition  of  the 
transmitter  input  plus  channel  noise.  This  is  accomplished 
by  operation  of  the  power  amplifier  in  a  linear  region  or  with 
saturating  power  amplifiers  using  constant  envelope  modulation 
techniques.  For  linear  systems,  multipath  fading  can  be  charac¬ 
terized  by  a  transfer  function  of  the  channel:  H(f,*t).  This 
function  is  the  two  dimensional  random  process  in  frequency  f 
and  time  t  that  is  observed  as  carrier  modulation  at  the  out¬ 
put  of  the  channel  when  sine  wave  excitation  at  the  carrier 
frequency  is  applied  to  the  channel  input.  For  any  continuous 
random  process,  we  can  determine  the  minimum  separation  required 
to  guarantee  decorrelation  with  respect  to  each  argument. 
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For  the  time  varying  transfer  function  H(f,*t),  let  tj  and  f d 
be  the  .decorrelation  separations  in  the  time  and  frequency 
variables,  respectively,  if  t<j  is  a  measure  of  the  time  de- 
correlation  in  seconds,  then 

°t Hz 

is  a  measure  of  the  fading  rate  or  bandwidth  of  the  random 
channel.  The  quantity  ot  is  often  referred  to  as  the  Doppler 
spread  because  it  is  a  measure  of  the  width  of  the  received 
spectrum  when  a  single  sine  wave  is  transmitted  through  the 
channel.  The  dual  relationship  for  the  frequency  decorrela¬ 
tion  fd  in  Hz  suggests  that  a  delay  variable 

<jf  *  seconds 

d 

defines  the  extent  of  the  multipath  delay.  The  quantity  af 
is  often  referred  to  as  the  multipath  delay  spread  as  it  is  a 
measure  of  the  width  of  the  received  process  in  the  time  domain 
when  a  single  impulse  function  is  transmitted  through  the 
channel. 

Typical  values  of  these  spread  factors  for  troposcatter 
communication  are 


TROPOSCATTER 

afc  *  1  Hz:  Doppler  Spread 

of  ~  10”7  seconds:  Delay  Spread 

where  the  symbol  ”  denotes  "on  the  order  of". 


The  spreads  can  be  defined  precisely  as  moments  of  spectra 
in  a  channel  model  which  assumes  Wide  Sense  Stationarity 


2-6 


• 

(WSS)  ip  the  time  variable  and  Uncorrelated  Scattering  (US) 
in  a  multipath  delay  variable.  This  WSSUS  model  and  the 
assumption  of  Gaussian  statistics  for  H(f,*t)  provide  a 
statistical  description  in  terms  of  a  single  two-dimensional 
correlation  function  of  the  random  process  H(f,  *t). 

This  characterization  has  been  quite  useful  and  accurate 
for  a  variety  of  radio  link  applications.  However  the  station- 
arity  and  Gaussian  assumptions  are  not  necessary  for  the 
utilization  of  adaptive  signal  processing  techniques  on  these - 
channels.  What  is  necessary  is  first  that  sufficient  time 
exist  to  "learn"  the  channel  characteristics  before  they  change 
and  second,  that  decorrelated  portions  of  the  frequency  band 
be  excited  such  that  a  diversity  effect  can  be  realized. 

These  conditions  are  reflected  in  the  following  two  relation¬ 
ships  in  terms  of  the  data  rate  R,  the  bandwidth  B,  and  two 
previously  defined  channel  factors,  VIZ; 

R(b/s)  >>  ot(Hz)  Learning  Constraint 

•  •••  ....  ^ 

*  *•*••••  •  •• 

B(Hz)  £  fd(Hz)  Diversity  Constraint. 

where  the  symbol  £  denotes  "  on  the  order  of  or  greater  than  " . 

The  learning  constraint  insures  that  sufficient  signal- 
to-noise  ratio  (SNR)  exists  for  reliable  communication  at 
rate  R  over  the  channel.  Clearly  if  R  -  afc,  the  channel 
would  change  before  significant  energy  for  measurement 
purposes  could  be  collected.  When  R  >>  o^,  the  received 
data  symbols  can  be  viewed  as  the  result  of  a  channel  sound- 
S  ing  signal  and  appropriate  processing  cam  generate  estimates 

of  the  channel  character  during  that  particular  stationary 
epoch.  The  signal  processing  techniques  in  an  adaptive  re¬ 
ceiver  do  not  necessarily  need  to  measure  the  channel  directly 
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in  the  optimization  of  the  receiver  but  the  requirements  on 
learning  are  approximately  the  same.  If  only  information  symbols 
are  used  in  the  sounding  signal,  the  learning  mode  is  referred 
to  as  Decision-Directed .  When  digital  symbols  known  to  both 
the  transmitter  and  receiver  are  employed,  the  learning  mode 
is  called  Reference-Directed .  An  important  advantage  of  digital 
systems  is  that  in  many  adaptive  communications  applications, 
adaptation  of  the  receiver  with  no  wasted  power  for  sounding 
signals  can  be  accomplished  using  the  decision-directed  mode. 

This  is  possible  in  digital  systems  because  of  the  finite  number 
of  parameters  or  levels  in  the  transmitted  source  symbols  and 
the  high  likelihood  that  receiver  decisions  are  correct. 

Diversity  in  fading  applications  is  used  to  provide  re¬ 
dundant  communications  channels  so  that  when  some  of  the  channels 
fade,  communication  will  still  be  possible  over  the  others  that 
are  not  in  a  fade.  Some  of  the  forms  of  diversity  employed 
are  space  using  multiple  antennas,  angle  of  arrival  using  mul¬ 
tiple  feedhorns,  polarization,  frequency,  and  time.  These 
diversity  techniques  are  sometimes  called  explicit  diversity 
because  of  their  externally  visible  nature.  An  alternate  form 
of  diversity  is  termed  implicit  diversity  because  the  channel 
itself  provides  redundancy.  In  order  to  capitalize  on  this 
implicit  diversity  for  added  protection,  proper  receiver  tech¬ 
niques  have  to  be  employed  to  correctly  assess  and  combine 
the  redundant  information.  The  potential  for  implicit  fre¬ 
quency  diversity  arises  because  different  parts  of  the  fre¬ 
quency  band  fade  independently.  Thus,  while  one  section  of 
the  band  may  be  in  a  deep  fade  the  remainder  can  be  used  for 
reliable  communications.  However,  if  the  transmitted  band¬ 
width  B  is  small  compared  to  the  frequency  decorrelation  in¬ 
terval  f^,  the  entire  band  will  fade  and  no  implicit  diversity 
can  result.  Thus,  the  second  constraint  B  £  f^  must  be  met  if 
an  implicity  diversity  gain  is  to  be  realized.  In  diversity 
systems  a  little  decorrelation  between  alternate  signal  paths 
can  provide  significant  diversity  gain. 
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It  is  easy  to  establish  that  for  the  transmission  rates 
of  interest,  troposcatter  systems  satisfy  both  the  learning 
and  diversity  constraint.  Digital  troposcatter  systems  are 
being  planned  with  data  rates  on  the  order  of  5  Mb/s  with 
approximately  the  same  bandwidth.  The  constraints  are  then 

R  *  5  x  10®  >>  ~  Is  Learning  Constraint 

B  *  5  x  10®  >  fd  ~  1.6  x  10®:  Diversity  Constraint 

This  means  that  adaptive  signal  processing  techniques 
can  be  applied  to  the  digital  troposcatter  receiver  with  the 
result  that  channel  variations  can  be  tracked  and  additional 
diversity  due  to  multipath  will  be  realized.  Note  this  result 
is  opposite  to  present  analog  technology  where  channel  tracking 
using  the  received  information  signal  is  not  possible  and  multi- 
path  degrades  rather  than  helps  performance.  This  concept  has 
been  successfully  exploited  to  realize  an  increase  in  digital 
rate  capability  of  about  an  order  of  magnitude.  A  digital  modem 
technique  1  *  '  ’  1  using  an  adaptive  decision- feedback 

equalizer  has  been  successfully  field  tested  at  digital  rates 
up  to  12.6  Mb/s  in  a  15  MHz  RF  bandwidth  allocation.  He  shall 
see  subsequently  that  not  only  does  this  new  technology  make 
angle  diversity  more  feasible  to  implement  but  angle  diversity 
provides  a  large  potential  for  bandwidth  conservation  in  the 
new  digital  systems. 

(2)  Conversion  to  Digital  Troposcatter  Systems:  Bandwidth 
Considerations 

The  majority  of  troposcatter  link  traffic  consists  of  speech 
transmissions.  In  an  analog  system  these  signals  are  frequency 
division  multiplexed  in  4  KHz  subchannels  to  produce  a  base¬ 
band  spectrum  of  approximately  60  KHz  +  N  x  4  KHz,  where  N  is 
the  number  of  voice  channels.  The  FM  system  expands  this  band¬ 
width  depending  on  the  selected  frequency  deviation.  If  we 
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let  Af  be  the  rms  frequency  deviation  per  channel,  the  peak 

frequency  deviation  is  /2n  Af  and  using  Carson’s  rule,  the 

rms 

FM  RF  bandwidth  is  approximately 


B__  ■  2  (4N+60+/2N  Af  )  KHz. 

rms 

For  a  120  channel  system  with  a  typical  150  KHz  rms 
frequency  deviation,  the  RF  bandwidth  is  5.7  MHz. 

A  120  channel  digital  system  with  64  Kb/s  Pulse  Code  Modu¬ 
lation  (PCM)  requires  a  data  rate  of  approximately  8  Mb/s  when 
framing  and  overhead  are  added.  This  example  illustrates  the  gen¬ 
eral  result  that  conversion  of  troposcatter  systems  from  analog  to 
digital  requires  modem  technology  with  better  than  1  bit/second/Hz 
packing  if  the  spectral  requirements  for  the  same  number  of 

message  channels  is  to  be  unchanged.  Present  modem  tech- 

f2  11  2  121 

nology1  *  '  *  achieves  1  bit/second/Hz  packing  in  the  tropo¬ 

scatter  channel  but  degrades  with  larger  densities.  At  1.3 
bit/second/Hz  the  degradation  is  about  1  dfi.  In  many  cases  the 
best  solution  is  to  combine  spectrum  efficient  modem  technology 
with  a  diversity  configuration  that  minimizes  system  spectral 
occupancy.  Angle  diversity  in  conjunction  with  space  diversity 
achieves  minimum  system  spectral  occupancy. 

(3)  Angle  Diversity  for  Bandwidth  Conservation 

Conventional  troposcatter  systems  use  two  antennas  per 
terminal  and  two  bandwidth  allocations  per  direction  to  achieve 
dual  space/dual  frequency  (2S/2F)  diversity.  Angle  diversity 
provides  redundant  paths  by  collecting  different  angle  of  arriv¬ 
al  signals  at  the  receiver.  A  dual  angle  diversity  system,  for 
example,  would  employ  a  two  horn  primary  feed  at  the  parabolic 
reflector  focus  instead  of  the  conventional  one  horn  primary 
feed.  In  the  frequency  diversity  system,  a  power  amplifier  is 
required  at  each  of  the  two  RF  carriers.  To  convert  from  dual 
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frequency  to  dual  angle ,  the  second  power  amplifier  is  tuned  to 
the  same  frequency  as  the  first  and  a  new  feed  structure  is 
provided.  The  resulting  dual  space/dual  angle  (2S/2A)  system 
now  requires  only  one  bandwidth  allocation  in  both  directions. 
Conversion  to  angle  diversity  halves  the  frequency  requirements 
of  present  2S/2F  quadruple  diversity  systems.  Since  both  diver¬ 
sity  configurations  have  4  diversity  outputs,  the  systems  are 
identical  from  the  input  to  the  receivers  down  through  the  modem 
to  the  data  output.  Figures  2.1  and  2.2  illustrate  the  2S/2F  and 
2S/2A  configurations.  The  2S/2A  system  does  not  require  a 
diplexer  as  .does  the  2S/2F  system,  but  an  additional  waveguide  run 
is  needed  from  the  feedhorn  to  the  receivers.  The  2S/2F  system 
requires  transmitters  tuned  to  two  frequencies,  a  single  horn  feed, 
and  diplexers  to  separate  the  frequencies.  The  2S/2A  system  has  both 
transmitters  tuned  to  the  same  frequency  and  a  dual  horn  feed.  The 
receivers  and  modem  are  identical  in  both  configurations.  The  second 
transmitter  in  the  2S/2A  configuratin  is  required  for  redundancy  in 
case  of  failure.  Its  on-the-air  use  provides  3  dB  more  scattered 
power  per  diversity  which  tends  to  offset  the  squint  loss  from 
the  elevated  beam  in  the  dual  angle  diversity  structure. 

This  savings  in  bandwidth  through  use  of  angle  diversity 
is  potentially  present  with  analog  troposcatter  systems,  but  it 
is  largely  dissipated  by  the  more  difficult  implementation  dis¬ 
advantages  of  analog  angle  diversity.  In  an  analog  system,  the 
time  delay  variation  between  the  main  and  elevated  beam  and  the 
additional  multipath  spread  of  the  elevated  beam  are  serious 
problems.  The  adaptive  processor  confeept  previously  described 
for  the  troposcatter  application  compensates  for  these  effects 
in  a  digital  system.  The  compensation  takes  the  form  of  adaptive 
equalization  using  the  received  signal  in  a  decision-directed 
technique . 
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DUAL  SPACE/DUAL  ANGLE  (2S/2A) 
QUADRUPLE  DIVERSITY  TROPOSCATTER  TERMINAL 


BASEBAND  ORDER  WIRE 


b.  Angle  Diversity  Desi 


The  short  term  (~  seconds)  fading  characteristics  of  tropo- 
scatter  signal  reception  are  mitigated  by  the  provision  of  redun¬ 
dant,  i.e.,  diversity,  channels  and  an  associated  combining  scheme 
to  utilize  the  stronger  received  signals  and  eliminate  the  weaker 
ones.  One  method  of  realizing  additional  diversity  channels  is 
to  use  multiple  feedhorns  at  the  focal  point  of  a  parabolic  re¬ 
flector  to  realize  multiple  angle  of  arrival  signals.  This  method 
•  •  •  •  •»■••.  •  •  «  •  •  •  . 

is  normally  only  employed  at  the  receiver  since  the  production  of 
multiple  transmit  beams  requires  either  additional  power  amplifiers 
or  a  reduction  in  power  with  power  splitters.  Thus,  conventional 
angle  diversity  has  one  transmit  beam  and  two  or  more  angle  of  ar¬ 
rival  receive  beams. 

One  of  the  first  questions  a  system  designer  may  ask  concerns 
the  relative  advantages  of  beams  spread  either  vertically  or  hor¬ 
izontally.  Since  the  diversity  advantage  stems  from  the  lack  of 
correlation  between  the  two  received  diversity  signals,  the  cor¬ 
relation  fall-off  as  a  function  of  squint  angle  in  the  vertical 
and  horizontal  directions  is  the  determining  factor  in  assessing 
the  performance  difference.  The  power  loss  fall-off  as  a  function 
of  squint  angle  is  of  the  same  order  for  the  two  directions  and 
therefore  is  not  significant  in  the  selection  process.  The  cor¬ 
relation  fall-off  as  a  function  of  squint  angle  can  be 
related  to  the  correlation  distance  in  the  plane  of  the  receive  * 

antenna.  The  correlation  distance  rc  may  be  defined  in  terms  of 
two  point  antennas,  one  located  on  the  great  circle  path  and  one 

at  a  distance  r  ,  where  the  received  processes  on  the  two  point 

c  % 

antennas  have  a  normalized  correlation  coefficient  of  1/e.  Clearly 

if  the  correlation  distance  is  much  greater  than  the  troposcatter 
system  parabolic  reflector  then  the  correlation  between  squinted 
beams  will  be  very  high  whereas  if  the  correlation  distance  is  • 
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• 

small  compared  to  the  parabolic  reflector  the  squinted  beam  cor¬ 
relation  will  be  small.  The  correlation  distance  is  defined  in 
terms  of  a  vector  u  in  the  plane  of  the  parabolic  reflector  as 

cp(u)  *  E(r)  E  (r  +  u) 

where  E(r)  is  the  received  electric  field  at  a  distance  r  from  the 
center  of  the  scattering  volume.  This  spatial  correlation  func- 
tion  is  related  to  the  scattering  volume  "size"  through  a'  three 
dimensional  Fourier  integral 

u.r 1 

PPf  ^  r  3 

cp(u)  *  C  JjJ  V(r* )  e  dr1. 

The  integration  is  over  the  scattering  volume  defined  by 

V(r’)  which  is  determined  by  the  antenna  gain  patterns  and  the 

refractive  index  spectrum  fall-off  as  a  function  of  scattering 

angle.  The  constant  k  is  equal  to  2nf  /c  where  f  the  carrier 

o  o 

frequency  and  c  is  the  speed  of  light.  Now  Fourier  variables  in 
one  domain  are  inversely  related  to  Fourier  variables  in  the 
other  domain,  e.g.,  a  short  time  pulse  corresponds  to  a  wide 
frequency  band.  Thus  the  correlation  distance  in  a  par¬ 
ticular  direction  at  the  receiver  is  inversely  related 
to  the  common  volume  size  in  that  same  direction.  In  Fig. 
2.3  we  examine  the  common  volume  dimensions  in  the  vertical  (side 
view)  and  horizontal  (top  view)  directions.  For  angle  diversity 
systems,  the  transmit  beamwidth  ft  is  typically  smaller  than  the 

*  minimum  scattering  angle  6  in  order  to  minimize  the  loss  associat- 

o 

ed  with  the  squinted  beam.  The  common  volume  "size"  in  the  hori¬ 
zontal  direction  is  then  limited  by  the  transmit  beamwidth  and  is 
(  on  the  order  of  Y^r  ■  ftr.  In  the  vertical  direction,  the  common 

volume  "size"  is  not  limited  by  the  beamwidth  as  much  as  the  re- 
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fractive  index  spectrum  fall-off  which  has  a  dependence  of  6 

where  0  is  the  scattering  angle.  Since  6  >  9q>  0  the  vertical 

common  volume  size  Y^r  is  larger  than  the  horizontal  common  volume 

size  YgT.  The  area  of  useful  scattering  energy  returned  to  the 

receiver  is  figuratively  shown  in  Fig.  2.3  as  a  shaded  area.  Since 

Y  >  Y_  it  follows  that  the  vertical  correlation  distance  will  be 
V  H 

less  than  the  horizontal  correlation  distance  for  narrow  beamwidth 
antennas.  Thus  the  vertical  angle  diversity  system  will  result  in 
a  lower  correlation  between  squinted  beams  than  a  horizontal  angle  * 
diversity  system. 

It  turns  out  that  as  the  transmit  beamwidth  is  increased 

such  that  fi  »  0  ,  the  situation  reverses  and  YTr  becomes  smaller 

than  Y  with  the  result  of  a  smaller  horizontal  correlation  dis- 
H 

tance  than  vertical.  This  result  is  consistent  with  **fat“  trans¬ 
mitter  beam  measurements  and  analysis  [2.4]  of  correlation  dis¬ 
tances  in  the  vertical  and  horizontal  planes.  However,  in  any 
practical  application,  the  transmit  beams  must  be  narrow  in  order 
to  provide  the  required  antenna  gain  for  successful  operation. 

Because  the  vertical  correlation  distance  is  smaller  than  the 
horizontal  correlation  distance  for  practical  systems,  the  angle 
diversity  design  should  utilize  a  vertical  squinted  beam  as  the 
first  additional  angle  of  arrival  signal.  The  use  of  more  angle 
of  arrival  signals  than  two  is  not  advantageous  both  because  of 
system  complexity  and  the  diminishing  return  from  additional 
diversity.  Thus  we  fix  the  system  design  as  a  dual  feedhorn  AD 
system  and  turn  our  attention  to  the  choice  of  vertical  squint 
angle  and  antenna  boresight/horizon  angle. 

The  angle  between  the  centerlines  of  the  antenna  patterns  of 
the  two  feedhorns  in  a  dual  AD  system  is  defined  as  the  squint 
angle.  When  the  squint  angle  is  appreciably  more  them  a  beamwidth. 


the  correlation  between  the  received  angle  diversity  signals  is 
small  but  the  relative  signal  loss  of  the  elevated  beam  may  be 
excessive  due  to  the  increased  scattering  angle.  Decreasing  the 
squint  angle  reduces  this  loss  but  increases  the  correlation  be¬ 
tween  diversity  signals.  The  effect  of  this  correlation  is  to 
reduce  the  diversity  gain.  It  becomes  necessary  to  distinguish 
between  narrowband  and  wideband  correlation  effects  on  diversity 
combining  loss.  When  the  diversity  signals  are  smeared  by  a 
"multipath  spread  that-  has-  width. on  the  .prder  of  the  reciprocal 
of  the  signal  bandwidth,  the  correlation  between  diversity  sig¬ 
nals  also  becomes  smeared.  Thus  the  diversity  combining  loss 
due  to  diversity  beam  correlation  is  less  for  wideband  (~  on 
the  order  of  a  reciprocal  multipath  spread)  than  for  narrowband 
(much  less  than  a  reciprocal  multipath  spread)  transmissions. 

In  the  propagation/modem  model  developed  under  this  program  the 
effects  of  correlation  on  wideband  signaling  is  correctly  taken 
into  account  through  the  introduction  of  a  cross  channel  multi- 
path  profile.  For  narrowband  signaling  the  signal  correlation 
loss  can  be  expressed  at  large  signal-to-noise  ratios  as 

2 

L  *  -5  LOG (1-p  )  (Narrowband  Correlation  Loss) 

where  p  is  the  normalized  correlation  coefficient  between  two 
complex  Gaussian  processes. 

A  simplified  model  developing  the  wideband  correlation  loss 
has  not  been  derived.  Instead  we  use  an  exact  analysis  requiring 
numerical  integration  over  the  common  volume.  The  propagation  pre¬ 
diction  model  and  the  modem  performance  model  described  subsequently 
in  Sections  3  and  4,  respectively  are  used  to  generate  modem  per¬ 
formance  curves  as  a  function  of  vertical  squint  angle.  The  perfor¬ 
mance  measure  used  is  outage  rate  defined  as  the  probability  that 
the  short  term  0.2  to  5  second  bit  error  rate  (BER)  exceeds  a  value  of 

10  .  The  BER  criterion  of  10~  has  been  selected  by  the  Defense 

Communications  Agency  because  at  this  level  PCM  voice  begins  to  de¬ 
grade.  The  use  of  outage  rate  has  been  established  as  a  more  mean¬ 
ingful  performance  measure  with  respect  to  digitized  voice  than 


2-18 


the  use  of  average  error  rate.  The  results  of  the  outage  rate 
calculations  for  an  example  C-band  and  L-band  (255-985  MHz)  link  are 
shown  in  Figures  2.4  and  2.5.  These  curves  reflect  the  effects  of  both 
squint  loss  and  diversity  correlation  as  a  function  of  squint 
angle  measured  in  beamwidths  (BW) .  The  two  examples  have  ap¬ 
proximately  the  same  main  beam  multipath  spread.  (20^7^0.45) . 
Figure  2.4  represents  the  RADC  Youngs town- Verona  test  link  while 
Figure  2.5  is  the  DCS  SAHIN  Tepesi-Yamanlar  link. 

The  wideband  correlation  loss  can  be  isolated  by  noting 
that  the  diversity  loss  in  the  absence  of  correlation  is  approx¬ 
imately  one-half  the  squint  loss. 

Performance  measures  with  unequal  branch  SNR  in  dB  are  well 
approximated  by  the  geometric  mean  of  the  SNR  or  the  arithmetic 
mean  in  dB.  Thus  by  comparing  the  diversity  loss  of  specific 
squint  angle  cases  with  a  case  with  zero  correlation  coefficient, 
the  wideband  correlation  loss  can  be  computed  from  the  curves  of 
Figures  2.4  and  2.5.  These  computations  are  summarized  in  Tables 
2.1  and  2.2.  Note  that  the  optimum  squint  angle  is  1  beamwidth 
for  narrowband  signaling  but  at  6.3  Mb/s  on  these  links  the  best 
performance  occurs  at  3/4  beamwidth  due  to  smaller  correlation 
effects  in  wideband  transmission.  Plots  of  the  correlation  coef¬ 
ficients  and  the  wideband  and  narrowband  correlation  losses  are 
given  in  Figures  2.6  and  2.7  for  the  two  examples .  The  reduced 
effect  of  correlation  between  angle  diversity  branches  with  wide¬ 
band  signaling  is  due  to  the  adaptive  equalizer  acting  to  combine 
individual  multipath  components  on  diversity  pairs.  These  in¬ 
dividual  multipath  components  have  less  correlation  than  the  cor¬ 
relation  from  the  cumulative  multipath  structure. 

In  these  examples  and  all  others  computed,  the  optimum 
squint  angle  was  determined  to  be  approximately  1  beamwidth  for 
narrowband  signaling  and  3/4  beamwidth  for  wideband  signaling. 

Thus,  it  is  clear  that  the  best  squint  angle  for  megabit  data  rates 
in  angle  diversity  systems  is  less  than  one  beamwidth.  Feedhorn 
design  constraints,  however,  make  it  difficult  to  achieve  squint 
angles  this  small.  Fortunately,  the  dB  loss  is  small  for  squint 
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Fig.  2.4  Outage  Rate  for  C-band  Link 
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Fig.  2.5  Outage  Rate  for  L-band  Link 
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angles  larger  than  3/4  beamwidth.  Because  of  the  feedhorn  design 
problem ,  this  result  has  important  practical  implications.  The 
feedhorn  design  for  the  RADC  test  link  resulted  in  the  measured 
squint  angles  given  in  Table  2.3.  Vertical  polarization  is  used 
for  signal  reception  in  the  RADC  tests.  The  additional  loss  for 
this  link  resulting  from  a  squint  angle  larger  than  3/4  beamwidth 
can  be  determined  from  Fig.  2.4  to  vary  between  0.3  and  0.8  dB. 


Table  2.3 

Test  Link  Squint  Angles 


Right  Antenna 

Left  Antenna 

H  Pol. 

V  Pol. 

H  Pol. 

V  Pol. 

Squint 

Angle 

0.65° 

0.73° 

0.55° 

0.68° 

Half-Power 

Beamwidth 

0.55° 

0.53° 

0.50° 

0.50° 

The  next  question  to  be  addressed  is  the  choice  of  antenna 
pointing  angle ,  i . e . ,  the  angle  between  the  main  beam  antenna 
pattern  centerline  (boresight)  and  the  horizon.  As  this  angle 
is  decreased,  the  relative  loss  of  the  elevated  beam  is  reduced 
but  more  of  the  main  beam  pattern  is  blocked  by  the  radio  horizon. 
Thus  an  optimum  pointing  angle  exists.  Since  operational  angle 

diversity  system  are  utilized  in  both  directions,  the  boresight/ 
horizon  angles  at  transmitter  and  receiver  are  constrained  to  be 
equal  in  the  optimization  search.  Using  a  squint  angle  of  one 
beamwidth  ,  the  main  beam  loss  for  a  conventional  system  was 
determined  as  a  function  of  antenna  pointing  angle  and  in  addi¬ 
tion  the  diversity  combining  loss  due  to  squint  loss  and  narrow¬ 
band  correlation  was  computed.  The  minimum  of  the  sum  of  these 

n - 

The  design  goal  for  the  feedhorn 
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losses  establishes  the  optimum  pointing  angles.  The  results 
for  3  example  systems  are  shown  in  Figs.  2.8  -  2.10.  In  these 
and  other  examples  the  optimum  bore sight/horizon  angle  for  the 
transmit  and  receiver  antennas  was  found  to  fall  between  1/4 
and  1/2  beamwidth  elevation  above  the  local  horizon.  This 
result  was  only  slightly  influenced  by  the  frequency  of  oper¬ 
ation  since  lower  frequency  systems  usually  use  larger  antennas. 

Note  that  the  major  effect  in  the  optimization  is  the  rapid 
increase  in  main  beam  loss  when  earth  blockage  begins  to  occur. 
For  this  reason  correlation  and  squint  loss  have  a  small  effect 
on  the  choice  of  optimum  pointing  angle.  This  large  dependence 
on  main  beam  blockage  also  establishes  that  the  optimum  pointing 
angles  for  conventional  and  angle  diversity  systems  are  approx¬ 
imately  the  same.  Although  other  researchers  have  suggested 
in  the  past  that  an  attempt  to  equalize  the  main  and  elevated 
beams  is  desirable,  these  results  show  that  the  main  beam  block¬ 
age  from  such  a  strategy  would  lead  to  an  overall  degradation. 
Under  homogeneous  atmospheric  conditions,  when  the  main  beam 
signal  is  maximized,  the  elevated  beam  signal  will  be  somewhat 
smaller.  The  frequent  occurrence  of  nonhomogeneous  atmospheric 
conditions  causes  the  elevated  beam  signal  to  be  stronger  than 
the  main  beam  signal.  This  phenomenon  makes  it  necessary  to 
aim  the  antenna  system  with  a  series  of  tests  extending  over  at 
least  several  days  in  order  to  preclude  pointing  the  antenna 
into  the  ground  as  a  result  of  nonhomogeneous  atmospheric  con¬ 
ditions.  It  is  likely  that  most  DCS  troposcatter  antenna  systems 
are  aimed  too  low  and  2  or  3  dB  could  be  gained  by  careful  aiming 
techniques . 


I 
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Antenna  Pointing  Angle  Optimization,  RADC  Test  Link 
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Fig.  2.9  Antenna  Pointing  Angle  Optimization,  S.  Tepesi-Yamanlar 


205  MILES  900  MHZ.  GG'AMTSNNAS 
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Pointing  Angle  Optimization,  Oslo-Kristiansand 


A  final  consideration  in  choice  of  antenna  pointing  angle 
is  the  effect  of  increased  multipath  in  both  beams  due  to  larger 
pointing  angles.  Because  digital  modems  can  utilize  the  multi- 
path  to  increase  the  effective  diversity  order  if  the  multipath 
spread  is  not  too  large,  it  is  advantageous  to  choose  a  somewhat 
larger  pointing  angle  if  the  nominal  multipath  spread  is  small 
compared  to  the  data  symbol  interval.  However,  for  L-band  systems 
the  multipath  spread,  particularly  in  the  elevated  beam,  may 
exceed  the  multipath  capabilities  of  the  digital  modem.  In  this 
situation  degraded  performance  due  to  intersymbol  interference 
results. 

The  results  from  the  numerical  integration  of  the  common 
volume  using  a  turbulent  scattering  hypothesis  are  shown  in 
Figs.  2.11  -  2.13  for  the  example  systems.  The  2a  multipath 
spread  of  the  product  beam  refers  to  the  cross  channel  multipath 
profile.  The  C-band  test  link  has  small  multipath  spread*  for 
data  rates  of  interest  and  because  of  the  small  dB  loss  asso- 
iated  with  larger  pointing  angles,  the  optimum  pointing  angle 
for  a  digital  AD  system  is  probably  very  close  to  the  optimum 
pointing  angle  for  a  conventional  2S/2F  analog  troposcatter 
system,  i.e.,  approximately  1/2  beamwidth.  For  the  L-band 
systems  there  is  considerably  more  multipath  and  a  steeper  slope 
with  pointing  angle.  Optimum  performance  of  a  digital  AD  sys¬ 
tem  at  these  frequencies  may  result  from  pointing  angles  some¬ 
what  less  than  1/4  beamwidth. 

In  summary  an  optimized  angle  diversity  system  utilizes 
a  squint  angle  of  approximately  3/4  beamwidth  with  a  boresight/ 
horizon  angle  approximately  equal  to  1/2  beamwidth  for  C-band 
digital  systems  and  somewhat  less  than  1/4  beamwidth  for  L-band 
systems.  Under  homogeneous  atmospheric  conditions  this  con¬ 
figuration  results  in  less  received  signal  power  in  the  el¬ 
evated  beam.  The  use  of  pointing  angles  at  much  less  than  1/4 
beamwidth  is  suboptimum  because  too  much  power  is  blocked  in  the 
main  beam  by  the  radio  horizon. 

Measured  results  from  C-band  systems  tend  .to  be  larqer  due  to 
the  inclusion  of  the  RF  terminal  equipment  filter  characteristics 
and  measurement  noise. 
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SECTION  3 


TROPOSCATTER  PATH  PREDICTION  TECHNIQUE 
a.  Introduction 

In  this  section  we  present  new  formulas  for  the  path-loss  and 
receiver  correlation  distances  for  several  values  of  the  spectrum 
slope  m.  For  m  »  11/3,  our  pathloss  formula  depends  on  certain 
standard  atmospheric  parameters,  while  for  m  *  5  it  corresponds 
directly  to  the  NBS  model.  This  work  forms  the  basis  for  the 
SIGNATRON  propagation  prediction  model.  This  model  and  its  as¬ 
sumptions  are  presented  here. 

The  necessary  background  in  the  properties  of  tropospheric 
transmission  is  provided  in  subsection  3.2.  Subsection  establi^ies  the 
basic  pathloss  formula  with  widebeam  antennas  and  discusses  the 
concept  of  aperture-to-medium  coupling  loss.  The  results  for 
narrow  beam  antennas  are  treated  in  subsection  3.4,  while  the 
path  loss  with  intermediate  beam-sizes  is  derived  in  subsection 
3.5.  The  section  concludes  with  a  summary  of  the  troposcatter 
prediction  model.  We  begin  with'  a  general  'discussion  of  previous 
work. 

Theory  and  practical  techniques  for  troposcatter  communica¬ 
tion  started  developing  around  1950  [Booker  and  Gordon,  1950]. 

Two  main  theories  have  been  proposed,  the  layer  reflection  theory  [Friis, 
et  al,  19573  and  the  turbulent  scattering  theory  based  on  the  turbulence 
results  developed  by  Obukhov  and  Kolmogorov  in  1941  [see  Tatarskii,  1971], 
Initially,  various  layer  reflection  theories  received  the  most  attention. 
Partly  because  the  predicted  behaviour  was  close  to  the  observed  linear 
dependence  of  the  scattering  cross-section  with  wavelength.  The  layer 
reflection  theory  also  predicts  a  dependence  on  the  scattering  angle 
close  to  the  empirically  determined  dependence  of  the  form  ?  . 

In  later  years,  experiments  at  higher  frequencies  and  with  narrower 
antenna  beams  [Eklund  and  Wicker ts,  1968;  Hardy  and  Katz  ,1969] 


TW 


have  indicated  the  applicability  of  the  turbulence  theory,  with  a 

-11/3  — *  •• 

scattering  angle  dependence  of  the  form  9  and  a  wavelength 

dependence  for  the  scattering  cross-section  of  the  form  X  ^3. 

In  practice  layer  reflection  and  scattering  will  exist  simulta¬ 
neously  with  layer  reflection  being  more  frequent  when  widebeam 
antennas  and  low  frequencies  are  used.  Experimental  results 
[Eklund  and  Wickerts  #1968]  indicate  that  the  wavelength  dependence 
can  vary  between  X  *  and  X3.  The  scattering  angle  dependence  has 
been  measured  to  lie  in  the  range  of  6  to  9  [Gjessing,  1969] . 

The  exponent  in  the  scattering  angle  is  frequently  called  the  re¬ 
fractive  index  spectrum  slope  since  it  also  indicates  the  fall-off 
of  the  wave  number  spectrum  of  the  refractive  index,  provided  the 
wavelength  is  in  the  so-called  inertial  subrange,  i.e. ,  between 
the  inner  and  outer  scale  of  turbulence. 

A  number  of  theoretical  and  computational  models  have  been 

developed  rsooker  and  Gordon, 1950 r Tatarskii, 1971?  Rice  et  al,1965; 

«  ...»*•«  ••  •  • 

Hartman  and  Wilker son, 1959;  Yeh,  1962].  Further  references  and  a 

discussion  of  these  models  can  be  found  in  Larsen, [1968]  or  Panter, 
[1972].  The  most  generally  accepted  model  is  that  of  NBS  [Rice, 
et.al,1965;  Hartman  and  Wilkerson,  1959]  which,-  based  on  a  large 
number  of  empirical  results,  assumes  the  refractive  index  spectrum 
slope  is  five.  Several  arguments  can  be  forwarded  in  support  of 
using  a  spectrum  slope  of  11/3  instead.  First,  the  trend  is  toward 
troposcatter  links  at  higher  frequencies  where  layer-reflection 
is  less  important.  Second,  on  links  with  both  turbulent  scat¬ 
ter  and  layer  reflection  the  turbulent  scatter  is  the  more  critical 
effect  since  during  the  worst  days,  communication  will  be  by  scat¬ 
ter  only  and  layer  reflection,  when  present,  tend  to  increase  the 
received  power. 
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Tropospheric  Scatter  Transmission 


Radiowave  propagation  in  'the  troposphere  is  determined  by 
the  refractive  index.  The  refractive  index  can  be  described  by 
its  short-term  mean  and  variance,  both  of  which  are  highly  depen¬ 
dent  on  geographic  location,  season,  and  time  of  day.  The  mean  re^ 
tractive  index  is  usually  linearly  decreasing  as  a  function  of 
height,  and  can  be  handled  by  defining  an  effective  earth  radius. 
This  well  known  technique  will  also  be  used  here,  but  it  is 
emphasized  that  for  situations  where  the  refractive  index  profile 
is  not  linear  the  effective  earth  radius  used  must  be  such  that 
the  correct  scattering  angle  is  found. 

The  refractive  index  flunctuations  are  characterized  by  the 
covariance  function.  If  n(r)  is  the  refractive  index  at  a  point 
with  coordinates  r  then  define 

n^(r)  =  n(r)  -  E[n(r)], 
and  its  covariance  function 

—  * 

■  ECv£i)  w]  • 


In  a  homogeneous  medium  this  covariance  is  only  a  function  of 
r^-  r2*  and  in  such  a  case  the  wavenumber  spectrum  is  defined 
*>Y 


♦a® 


(2tt) 


i  err  ,  *  ik*r  .3 

r  jj  CD  (r)  e - d  r. 
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Often  the  structure  function  is  defined  instead  of  the  covariance 
function  [Tatarskii,  1971] , however  the  present  discussion  will  be 
simplified  by  assuming  that  the  volume  of  the  troposphere  concerned  can 
be  considered  homogeneous  and  isotropic.  For  the  isotropic  atmosphere 
the  wavenumber  spectrum  is  only  a  function  of  the  magnitude  k  of 


the  wavenumber  vector  k.  The  wavenumber  spectrum  is  important 
since  it  has  been  shown  (Tatarskii ,  1971]  that  the  scattering  cross- 
section  (a  )  of  a  small  volume  (dV)  is 

S 

a  =  8 it2  k4  $  (2k  sin  |)  dV, 

s  n  z 

where  8  is  the  scattering  angle.  Polarization  losses  are 
usually  treated  separately  and  have  been  neglected  in  this 
equation.  The  wavenumber  spectrum  is  usually  modeled  by  the 
von  Karman  spectrum  [Tatarskii ,1971] 


§  OO 
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(i 


m/2 


k2r 


2\ 

) 


m  >  3 
(3.1) 


1£  (k)  is  plotted  on  a  doubly  logarithmic  scale  against  the  wave- 

number  k,  the  slope  falls  off  as  m,so  that  m  is  called  the  spectrum 
slope,  m  has  to  be  greater  than  3  for  *n  (k)  to  be  integrable  over  the 
three  dimensional  wavenumber  space.  .Values  of  m  smaller  than  3  are 
possible  only  in  a  limited  wavelength  region.  The  variance  of  the  re¬ 
fractive  index  is  then  limited  only  by  the  small  scale  structure 

of  the  turbulence.  It  will  be  assumed  here  that  the  small  scale 

2 

turbulence  can  be  ignored,  so  that  m  is  larger  than  three.  <rn 

is  the  variance  of  the  refractive  index  and  r  is  the  correla- 

o 

tion  distance  of  the  turbulence.  The  wavelength  \  »  2rr/k  is 
assumed  to  be  in  the  inertial  subrange: 


* 


l  «  X  «  L  ,  ft 

o  o 

where  and  are  the  inner  and  outer  scales  of  turbulence, 

respectively,  r  in  (1)  can  usually  be  identified  with  L  . 

o  o 
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The  power  received  from  a  small  scatterer  at  a  point  r  and  with 


the  cross-section  a  is 

s 


dP(r )  =  P  • 
R  T 


G(r) 
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a  (r) 


4ttRt2  (r )  4ttRr2  (r ) 


4tt 


where 


GT(r)  GR(r>) 


rt(-)(rr(-)) 


=  Transmitted  power  (received  power) 

*  Gain  of  transmitter  (receiver) antenna  toward 
the  scattering  point  r. 

*  Distance  from  transmitter  (receiver) 
to  the  scattering  point  r. 


Using  the  egression  for  the  scattering  cross-section  and  the 


j  / 


von  Kannan  spectrum  (1)  assuming 


kr  »  1 
o 


it  is  found  that  the  received  power  is 
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where 


(3.2) 


c  -  ^  ro3-“  K2-"  r(|)/(^  r  (ajl))  .  <3.3» 

and  the  range  of  integration^,  is  the  total  common  volume.  It 

2 

is  assumed  that  <j  and  r  do  not  vary  appreciably  throughtout 

n  o 

the  common  volume. 

While  the  Kolmogorov-Obukhov  theory  of  turbulence  predicts 
m«ll/3,  the  above  expression  can  be  used  formally  fbr. other  values 
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of  m.  In  particular,  setting  np>5  yields  the  expression  used  by 

the  NBS  model  tRice.et  al.,1965].  In  order  to  predict  the  re- 

2 

ceived  power  P  it  is  necessary  to  specify  m,  r  and  cr  .  As  we 
k  on 

assume  r  equal  to  L  ,  we  can  use  the  value  of  L  found  by 
o  o  o 

Fried  [1967] .  He  found  the  following  approximate  expression  for 
the  outer  scale  of  turbulence. 


L  2  Vh  *1  meter 
o 

where  h  is  the  height  of  the  air  mass.  It  must  be  emphasized 

that  considerable  variation  from  this  value  of  Lq  can  be  expected 

under  different  atmospheric  conditions.  Measurements  and  predic- 
2 

tions  of  ern  have  been  presented  by  numerous  authors  [Hardy  and 

Katz, 1969  ;  Fried, 1967  ;  Gossard, 1977] .  At  radio  frequencies, the 
2 

variance  is  strongly  dependent  on  the  humidity.  For  the  dry  at¬ 
mosphere,  measurements  at  radio-  and  at  optical- frequencies  yield 
— *  ia.ntic.X  v,1Ues  o£  *  Since  ^ 

the  worst  case  for  radio  wave  scatter  measurements,  it  will  often 
be  possible  to  apply  optical  results  to  the  study  of  troposcatter 
at  radio  frequencies.  This  is  useful  since  many  more  measurements 
have  been  made  of  the  optical  refractive  index  than  of  the  radio 
refractive  index.  Fried  [1967]  suggested  the  following  model  of 
the  optical  refractive  index  variance: 


an,opt“  6‘7 


-14  -h/3200  m 

10  e 


New  results  indicate  a  stronger  height  dependence  [Gossard, 1977  ; 
Brookner  ,  1971] , particularly  for  continental  air  masses. 


2 

a  is  the  main  source  of  variability  in  received  signal 

n  2 

levels,  so  a  complete  description  of  cn  variation  versus  weather, 
time-of-day,  season,  location,  etc.,  is  therefore  desirable.  Un¬ 
fortunately,  only  sporadic  measurements  have  been  made.  Gossard, 
[1977]  has  calculated  both  optical  and  radio  refractive  index 
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variations  of  several  classes  of  air  masses.  Using  this  method, 

2  2 

profiles  of  a.  (or  the  structure  constant  C  )  cam  be  calculated 
n  n 

as  a  function  of  geographic  location,  season,  and  time-of-day. 

2 

Radiosonde  data  can  also  be  used  to  predict  the  variance  or*  for 

a  troposcatter  link  using  a  formula  by  Tatarskii,  [1971], 

2 

[Sirkis,  1971].  Typical  values  of  a  at  heights  of  2  km  or  below 

-14  -13  “ 

are  in  the  range  10  to  10 

c.  The  Basic  Pathloss 

The  approximate  pathloss  with  omnidirectional  antenna  is 

now  derived.  The  condition  required  for  the  validity  of  the 

approximation  is  that  the  minimum  scattering  angle  9  be  small. 

s 

The  geometry  of  the  path  is  shown  in  Fig.  3.1,  where  also  the 
necessary  notation  is  defined.  The  integral  over  the  common 
volume  is  evaluated  by  first  integrating  parallel  to  the  z-axis 
(perpendicular  to  the  great  circle  plane  containing  transmitter 
and  receiver ) . 
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All  subscript  o  parameters  are  in  the  great  circle  plane  in 
Fig.  3.1.  When  the  scattering  angle  is  small  the  following  ap¬ 
proximation  is  valid: 
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where 


RoRRoT 

RORWOT 


The  integral  is  now  reduced  to 

n/2  tt/2 

J  d«r  d%  [(RoR+Ror)*ln  9o(2  *in  t)*"1]'1  . 
“to  “He  (3.6) 


VPT  - c  Bd  nr) 


where  B(x,y)  is  the  Beta  function. 


v)  r <*>r (y) 

B(x'y)  r(x+y)  * 

Since  the  scattering  angle  is  small*  so  that  sin? (r)«0 (r) ,  and 
m  >  3,  we  find  after  some  simple  calculations  that 


C  B 


(1  m-l\ 
\2*  2  ) 


u  /o  ~  *  *  0 

“  d(m-l)  (m-2) 


*2— m 


2  3-m  .  2-m  . 

n  o  4(m-l)(m-2) 


im-3I 


(6g"m)/d 


The  parameters  in  this  equation  are 
2 


'n 


(3.7) 

(3.8) 


>  variance  of  refractive  index  fluctuations 

r  *  correlation  distance  of  refractive  index 
o 

k  »  2tt/x  °  wave  number 
m  “  spectrum  slope  (m  >  3) 

9a  »  aRG+am^“  minimum  scattering  angle  (angular  distance  in 


TG 


distance  between  the  terminals; 


radians) 


/ 

3-9 


J 


The  result  (8)  is  valid  for  a  small  scattering  angle,  a  condition, 

which  will  almost  always  be  satisfied  in  practice.  Height  varia- 
2 

tions  in  o  have  been  assumed  weak  enough  so  that  C  can  be  taken 
n 

outside  the  integral  sign  in  (3.6)  and  becomes  a  coefficient  in  equa' 
tion  (3.7).  Comparison  of  (3.8)  with  an  exact  integration  of  (3.4) 
the  special  case  where  m  is  an  even  integer  and  the  horizon  eleva¬ 
tion  angles  are  zero  [Waterman, 1958]  shows  that  the  approximations 
we  have  made  are  exceedingly  accurate  as  long  as  the  distance  d  is 
much  smaller  than  the  effective  earth  radius.  Since  this  condition 
is  generally  true,  our  result  can  be  used  for  all  practical  tropo- 
scatter  distances. 

The  initial  assumption  of  omnidirectional  antennas  can  be  re¬ 
laxed  to  generalize  (3.8)  to  antenna  beams  wide  enough  to  encompass 
the  main  scattering  volume,  except  that  antenna  gain  is  then  no 
longer  unity  and  has  to  be  taken  into  account. 

The  importance  of  the  spectrum  slope  is  evident  from  the  re* 
suit.  As  k  *  2 rrf / c ,  it  is  clear,  for  instance,  that  setting  m*5 
yields  the  cubic  frequency  dependence  assumed  in  most  models  of 
tropospheric  scatter.  Since  the  frequency  dependence  with  m  *  11/3 
is  much  weaker  (f  5/^3) ,  it  is  not  surprising  the  turbulent  scatter 
provides  the  dominant  transmission  at  higher  frequencies.  Figure 

3.2  illustrates  some  typical  curves  calculated  for  m  *  5  using  the 

short  distance  NBS  formula  [Rice,  et  al.,1965]  and  our  (3.8)  for 

2  -14 

m  «  11/3  using  a  5  •  10  and  r  »  L  *  70  m.  It  can  be 

n  o  o 

seen  that  at  high  frequencies  or  long  distance  turbulent  scat¬ 
ter  tends  to  be  more  important  than  layer  reflection  or  feuillet 
scattering.  This  leads  to  the  important  conclusion  that,  in 
general,  it  is  not  possible  to  use  a  single  model  (fixed  m)  to 
compare  the  performances  at  different  frequencies  or  distances. 
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Scatter  Theory  with  Typical  Parameters 


Another  conclusion  is  that,  if  only  one  model  were  to  be  used, 
it  is  better  to  base  a  troposcatter  link  design  on  the  turbulent 
scatter  model.  This  model  is  accurate  at  high  frequencies  and 
is  useful  at  lower  frequencies  to  estimate  the  performance  at 
the  times  when  atmospheric  layering  is  absent.  This  represents 
the  kind  of  "worst  case”  situation  that  can  be  expected  on  a 
dry  winter  day,  for  instance. 

If  a  beam  is  narrow  then  the  pathloss  calculated  above  is 
no  longer  valid  since  not  all  potential  scatters  are  illuminated 
by  the  beam.  It  is  customary  to  account  for  this  by  introducing 
the  aperture-to-medium  coupling  loss.  This  coupling  loss  is  de¬ 
fined  by 

Coupling  loss  in  dB  »  total  path  loss  in  dB 

+  antenna  gains  in  dB 

-  basic  pathloss  in  dB  with  omni¬ 
directional  antennas. 

It  is  important  to  realize  that  the  coupling  loss  and  the 
basic  pathloss  are  complementary  and  must  relate  to  the  total  path 
loss  for  one  specific  model.  It  is  a  common  mistake  to  “calculate" 
the  total  path  loss  from  the  coupling  loss  and  basic  path  loss  mixing 
values  from  different  tropospheric  models.  It  is  also  necessary 
to  take  care  when  measuring  the  coupling  loss  directly  by  com¬ 
paring  transmission  with  a  small  and  a  large  antenna  on  a  given 
link  since  the  result  of  such  measurements  depends  on  highly 
variable  atmospheric  conditions  outside  the  smaller  common  volume 
of  interest.  While  these  measurements  can  provide  useful  statis¬ 
tical  information,  it  is  not  possible  to  identify  the  measured 
results  directly  with  the  coupling  loss  of  a  pathloss  prediction 
model  when  the  assumption  of  homogeneity  may  not  be  valid. 

In  the  following,  approximate  formulas  will  be  presented 
for  the  coupling  loss  for  different  values  of  m. 
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cl.  Formulas  for  Extremely  Narrow  Antenna  Beams 

When  one  or  both  of  the  antenna  beam  widths  tend  to  zero, 
then  it  is  again  possible  to  obtain  simple  formulas.  Suppose 
first  that  the  transmitting  beam  is  wide  while  the  receiving 
beam  is  narrow.  The  assumption  of  a  narrow  beam  is  seldom 
satisfied  to  the  degree  required  here,  but  the  results  are  still 
useful.  They  provide  a  simple  bound  on  performance,  and  can 
be  compared  with  similar  results  by  other  authors  [6 jessing, 1969; 
waterman, 1958;  Gjessing  and  McCormack, 1974] .  The  derivations 
of  the  results  below  follow  the  same  lines  as  exemplified  by  the 
derivation  in  Section  III. 

we  first  evaluate  expressions  for  the  horizontal  and  verti¬ 
cal  correlation  distances  in  the  plane  of  the  receiver  aperture. 
These  quantities,  while  of  interest  in  themselves,  can  be  re¬ 
lated  to  the  coupling  loss.  Instead  of  viewing  the  coupling  loss 
as  caused  by  neglecting  to  illuminate  certain  scatterers,  it  is 
now  viewed  as  a  result  of  the  incoherency  of  the  received  field- 
over  the  aperture  A^.  Let  the  received  field  Eq (r ) ,  where  r  is 
a  point  in  the  aperture,  have  the  correlation  function  p ( r  -  r'), 

P (r,  r')  *  E^Eo (r)  Eo(r')}  . 

The  received  power  is  then  proportional  to 

pR  -  II  a2£  II  a2*'  p  <«■  >  • 

**  *R 

The  loss  due  to  decorrelation  of  the  field  over  the  aperture  is 
then 
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-  10  log  “ 


Ajp(0) 


,'rrr  d2r  d2r ’  p(r-r‘). 

il  M  il  11  ™  “  * 


This  can  be  shown  to  be  equivalent  to  the  coupling  loss  defined 
earlier,  but  for  the  receiving  aperture  only.  A  crude  approxima¬ 
tion  to  L  can  be  found  in  terms  of  the  horizontal  and  vertical 
c 

correlation  distances  and  defined  below, 

Lc  -  -  10  log  (Lg  Iy^) . 

The  horizontal  and  vertical  correlation  distances  are  now 
found.  Let 

PHU)  and  PVU) 

be  the  horizontal  and  vertical  correlation  functions  at  the  re¬ 
ceiver.  The  horizontal  correlation  distances  L  is  then  defined 

n 


lh«  J  Pa(*)dt/PH(0>. 


We  find  by  an  integration  analogous  to  the  one  leading  to  (3.8) , 


_X_  m-2  r  1  /aTGV 

h  *  8  _/l  m-l\  L  n+m-1  \  9  .  • 

8  BV2'"T";  n*0  8 


(3.9) 


Note  that  for  a  troposcatter  link  with  equal  take-off  angles  at 
the  transmitter  and  receiver  we  have 


«w/9.  -  1/2  ■ 
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The  evaluation  of  the  vertical  correlation  distance  should 
take  into  account  the  pointing  angle  of  the  antenna.  If  the  re¬ 
ceiver  beam  is  steered  for  maximum  power  a  simple  result  is  ob¬ 
tained: 


(m-2) . 


(3.10) 


The  coupling  loss  on  a  link  with  a  widebeam  transmitter  and  a 
narrow  beam  receiver  cam  be  expressed  in  the  form 


L  -  10  log  [-*  — ^ 
c  L4tt  m 


/i  a=i' 

Bn2*  2  ) 


d  a 
s  RG 


(3.11) 


where  G  is  the  (large)  gain  of  the  receiving  aperture  (natural 
K 

units) .  These  asymptotic  results  (in  the  sense  of  an  essentially 
infinitely  large  receiving  aperture)  may  be  compared  with  equivalent 
results  by  other  authors.  Gjessing  and  McCormick ,  [1974]  report  sim¬ 
ilar  expressions  under  the  additional  assumption  of  zero  degree  take¬ 
off  angles  so  that 


0  =  2  a  =  2  a  =  — 

S  TG  RG  R 

e 

(R  is  the  effective  earth  radius) .  Their  results  are  in  rea- 
© 

sonable  agreement  with  ours  (3.9),  (3.10) ,  and  (3.11)  (less  than  3  dB 

difference  in  L  ) . 

c 


If  the  transmitting  antenna  beam  is  also  narrow,  then  the 
parallel  results  to  the  equations  (3.9),  (3.10)  and  (3.11),  become 

H  ' 

\  -  J-  (m-1) 


and 


(m— 1) (m-2) 


(3.12) 


In  the  expression  (3.12)  for  the  cowling  loss  Lc  it  has  been 
assumed  that  BR>/eT<H  <  where  SRH  and  STH  are 

the  horizontal  beamwidths  of  receiver  and  transmitter .respec¬ 
tively.  If  this  is  not  the  case  the  role  of  receiver  and  trans¬ 
mitter  is  simply  reversed  before  (12)  is  used.  G  is  again 

K 

the  gain  of  the  receiving  aperture,  but  GT  v  is  the  gain  of  the 
transmitter  in  the  vertical  plane  only. 


e.  Approximate  Formulas  for  Intermediate  Beamwidths 


So  far  we  have  only  found  accurate  formulas  for  two  ex¬ 
treme  cases,  very  wide  beams,  and  very  narrow  beams.  For  the 
majority  of  troposcatter  links  encountered  in  practice  neither 
of  these  situations  occur.  It  is  therefore  important  to  tackle 
the  more  difficult  problem  of  finding  the  pathloss  for  beams 
of  a  more  realistic  size.  Exact  results  can  only  be  found  by 
a  numerical  integration,  but  useful  approximations  will  be  found 
on  the  basis  of  ideal  rectangular  beam  shapes.. 


Let  V'  8t  H'  Pr  v'  Pr  H  denote  the  vertical  and  horizontal 
beamwidths  of  the  transmitter  and  receiver  antennas.  For  ideal 
rectangular  beams  the  limits  of  integrals  in  (5)  are  changed  so 


that 


PD/P_  - 
R  T 
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dz  1(WZ) 

"zo(vV 


(3.13) 
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where 


I(Wz) 

is  the  integrand  in  (5)  and 


SW  *  2  min  At  eT,H-  hor  6r.h] 


-  2  si~n~;o(aR|aT)  min  <BT.H  sinVBT.R3inaT> 


(3.14) 


This  integral  is  now  approximated  in  some  special  cases. 


A.  Wide  Horizontal  Beamwidths 


In  this  case  it  can  be  assumed  that  z  ~  •,  and  it  is  then 

o 


found  that 


P  /p  r 


c  B  (j^)  .  2-m  -B-  -  B- 


R'  T  d(m-l)  (m-2) 


fl  2-»  p  ^ V\ 

s  Fm-2  \  A  *  fl_  • 


(3.15) 


where  we  have  defined  the  function 

FN(x,y)  -  1  -  (l+x)"N-(l+y)“N  +  (l+x+y)"N. 

(3.15)  clearly  reduces  to  (3.7)  for  very  wide  vertical  beamwidths. 
The  coupling  loss  is 


6  0 

„  /PT.V  PR. VN 

-10  109  P»-2  ~f)  ■ 
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Narrow  Horizontal  BeamwidthSi 


Rather  than  attempting  an  exact  integration  we  use  the 
simple  bounds 


srna  •  sma 
_ R  _ T 


+  sina_ 


sinaR 


_  2Zo  _  3inaR  '  *inaT 

d0  ,inCSi  +  SinaT 


Invoking  the  small  angle  approximation  and  the  assumption  that 
3  is  small  the  upper  bound  yields 

XI 


2f3„  0 


p  /p  SC.3  S-r-  p  (lit* 

PR/PT  d  m(m-l)  m-l  \  0 


R.V' 

0  i 


The  function  F  .  (x,y)  is  the  same  as  defined  above.  The  error 
m-l 

in  this  bound  is  less  than  3  dB ,  and  this  error  can  be  sub¬ 
stantially  reduced  by  using  an  approximation  between  the  lower 
and  upper  bounds.  The  correct  approximation  is  the  one  which  is 
equivalent  to  (3.12)  for  very  narrow  beams.  This  approximation  is 

K  A  Art  .1— m  n  A 


Since 


PT  2‘max(cLTG,aRG)  d  m(m-l)  m- 


+  “rg 
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s  s 


the  coupling  loss  is  found  to  be  approximately 


Lc  -  10  log 
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This  reduces  to  (12)  for  very  narrow  beams.  axG_aRG  ma^  in 
practice  be  approximated  by 


aTG**aRG 


5vV 


+  ^T0"®RG** 


where  h^  and  HR  are  the  altitudes  of  the  antennas,  and  @TG,  0RQ 
are  the  horizon  grazing  angles  relative  to  horizontal. 
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C.  Narrow  Horizontal  Beamwidths, 


?,H  >'“"^R,  H 


In  this  case  the  horizontal  dimension  is  essentially  deter¬ 
mined  by  the  receiver  beamwidth,  and  the  following  approxima¬ 
tion  is  used: 

2  Zo  RoR  PR,H* 

This  case  is  more  complicated,  and  we  just  state  the  result  of 
the  integration: 


P_/P_ 
R  T 


d(ra-l)  s 


0  k 
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+k-l\  0  / 

S  L. 


,  6r 
(1  +  ^ 


1-m-k 


6  k  /  B  1— m— k  B  +  B  1— m— k  \ 

T  \  s  s  /_ 

(3.17) 


The  corresponding  coupling  loss  also  reduces  to  (3.12)  for  narrow 
beams. 

D.  Combined  Formulas 

The  formulas  may  be  combined  to  yield  useful  approximations 
over  the  entire  range  of  parameter.  Consider  first  the  case 
where  B-,  „  »  P_  Figure  3.3  shows  the  approximations  in  A.  and 
C. above.  Also  shown  is  a  combined  curve  calculated  by  simply 
adding  the  coupling  losses  (in  natural  units) .  This  is  a  rather  arbi¬ 
trary  approximation  which ,  however ,  is  asymptotically  good .  For  the 
case  illustrated  in  Fig.  3.3  it  is  seen  that  the  approximation  is  excel¬ 
lent  when  compared  with  a  numerical  integration  over  the  common  volume. 
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Comparison  of  Derived  Coupling  Loss  with  Numerical  Integration  -  Asymmetric  Path 


This  integration  used  the  non-rectangular  beamshape  associated 
with  a  parabolic  antenna.  Figure  3.4  shows  a  similar  comparison 
for  the  case  of  antennas  of  equal  size.  It  is  seen  again  that 
a  simple  addition  of  the  coupling  losses  yield  an  excellent 
approximation  to  the  more  exact  numerical  integration.  In  Figs. 
3.3  and  3.4  the  beamwidths  used  in  the  formulas  and  the  size  of 
the  parabolic  antennas  are  related  such  that  the  correct  asymp¬ 
tote  is  achieved.  This  relation  of  beamwidth  B  to  antenna 
diameter  D  is  simply 


where  A  is  the  wavelength. 

In  the  turbulent  scatter  model  of  tropospheric  scatter 
communication  developed  here,  the  total  pathloss  on  such 
a  link  is  divided  into  two  terms,  the  basic  pathloss  valid  for 
wide  beams,  and  the  aperture-to-medium  coupling  loss  accounting 
for  the  additional  loss  of  narrow  beams.  It  is  important  that 
basic  pathloss  and  coupling  loss  relate  to  a  single  model  of 
the  scatter  mechanism  in  order  to  yield  a  meaningful  total  path- 
loss. 

We  have  shown  how  the  atmospheric  parameters  can  be 
related  to  the  turbulent  scatter  model.  Combined  with  the 
closed-formed  expressions  for  the  pathloss  derived  in  this 
paper,  this  results  in  a  useful  method  for  predicting  the  per¬ 
formance  of  a  link.  However,  more  work  is  required  to  relate 

2 

the  important  parameters  on  and  rg  to  measurable  or  predict¬ 
able  atmospheric  parameters.  The  method  of  Gossard  [1977]  may 
be  used  to  find  the  refractive  index  variance  for  many  different 
climates.  Radiosonde  data  may  also  be  used  to  determine  this 
variance,  which  can  then  be  used  in  our  model  to  predict  com¬ 
munication  performance.  We  point  out,  however,  that  formation 
of  atmospheric  layers  can  partly  invalidate  the  turbulent  scat- 


ter  assumption  used  here.  At  lower  frequencies  (<  1  GHz) ,  a 
separate  analysis  of  layer  effects  is  required  to  complete  a 
pathloss  prediction  program.  The  NBS  method  [Rice  et  al,  1965] 
accomplishes  this  by  using  empirical  data.  However,  the  turbu¬ 
lent  scatter  model  is  most  useful  at  higher  frequencies,  and 
even  at  lower  frequencies  when  concerned  primarily  with  the 
weak-signal  situation. 

The  asymptotic  form  of  the  coupling  loss  is  sometimes 
used  to  approximate  the  actual  loss  [Gjessing  and  McCormick, 
1974].  Our  results  indicate  that  the  coupling  loss  is  very 
slow  in  approaching  the  asymptote  and  that  use  of  the  asymptotic 
result  is  much  too  optimistic  (see  Figs.  3.3  and  3.4).  The 
formulas  derived  in  subsection  3.5  provide  an  alternative  to 
using  the  asymptotic  form  without  this  deficiency. 

f .  Summary  of  Troposcatter  Model 

In  this  subsection  we  summarize  the  troposcatter  model 
developed  by  SIGNATRON,  and  its  relation  to  atmospheric  and 
climatic  factors.  The  dependence  of  path  loss  on  key  para¬ 
meters  is  also  developed. 

(1)  Path  Geometry 

A  typical  path  is  shown  in  Fig.  3.5  with  some  of  the 
parameters  involved.  The  parameters  are 


•  Frequency  f 

•  Distance  d 

•  Heights  h„,  h_  of  transmitter  and  receiver 
antennas  above  mean  Sea  Level  (MSL) . 

•  Horizon  elevations  above  horizontal  (0tG»9Rg) 
at  transmitter  and  receiver.  These  angles 
can  depend  on  the  atmospheric  conditions.  - 

•  Antenna  boresight  elevations  above  horizontal 
(9t»@r)  at  transmitter  and  receiver. 

Note  that  9t-9tg  and  6R-6RG  are  then 
the  elevation  or  the  antennas  above  the 


horizon.  The  definition  of  "boresight"  need 
not  be  the  direction  of  maximum  gain,  but 
may  be  defined  from  mechanical  consideration 
(e.g. ,  the  normal  to  a  phased  array). 

e  Azimuth  pointing  angles  67,4it  if  horizontal 

diversity  is  employed,  or  beam  swinging  experi¬ 
ment  performed. 

e  Gain  patterns  Gf(6, 4),  Gr(0,4)  of  trans¬ 
mitting  and  receiving  antennas.  Determination  of. 
the  gain  patterns  requires  a  number  of  additional 
parameters  -  size  and  type  of  antennas,  height 
above  ground,  ground  profile  and  reflectivity. 

The  patterns  should  be  referenced  with  respect  to 
the  boresight  direction.  Both  amplitude  and 
.  phase  of  the  voltage  gain  patterns  are  required 
when  the  cross  correlation  between  two  (diver- 
sity)  paths  is  to  be  calculated. 

e  Polarization  of  transmitter  and  receiver. 


(2)  Turbulent  Scatter  Model 

It  is  well  known  that  the  power  received  from  a  small 


scattering 

volume 

can  be  written 

PR 

*  PT 

6T  as  X2<5r  .-1 
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where 

PT 

- 

• 

Transmitter  power 

gT 

« 

Power  Gain  of  transmitting  antenna 

*T 

M 

Distance  from  transmitter  to  scatterer 

*• 

• 

Scattering  cross-section  of  the  turbulent 
volume 

rr 

- 

Distance  from  scatterer  to  receiver 

1 

m 

Wavelength 

Gr 

- 

Power  gain  of  receiving  antenna 

Aa 

- 

Atmospheric  attenuation 
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Th«  above  transmission  equation  also  serves  as  a  definition  of 
the  scattering  cross-section  (some  authors  define  ae/4v  as 
the  scattering  cross-section).  Tatarski  (1971)  has  derived  the 
cross-section  of  turbulent  scatter  with  a  given  wave  number 
spectrum# 


•*  *  8,2  **  *n  (2*  sin  § )  sin2  x  dv#  *0  «1«  LQ. 


In  this  equation# 


k  *  2ir/X  is  the  wave  number# 

*n(k)  is  the  locally  homogeneous  and  isotropic 

wave  spectrum 

6  is  the  scattering  angle 

sin2x  accounts  for  loss  due  to  polarization  mis* 

match  (usually  negligible) 

X  angle  between  incoming  electric  field  and 

direction  of  propagation 

dV  is  the  infinitesimal  volume. 


The  equation  is  derived  under  the  assumption 


Lo  «  /*Rr  ,  /lR^ 


If  a  spectrum  of  the  von  Karman  type  is  assumed#  it  is 
found 


*.  ■  9n  8 


/^4-m  3-m  Lil 


2  sin  ^ 


r0  is  approximately  the  correlation  distance  L©  of  tne  tur¬ 
bulence.  Polarization  losses  have  oeen  ignored  in  this 
equation.  For  m  *  11/3  we  get  the  cross-section  prediction  by 


the  Kolmogorov-Obukhov  theory,  while  for  m  =  5  we  get  the  for¬ 
mula  on  which  the  NBS  method  [Rice,  et.  al,  1967)  is  based. 
Measurements  agree  with  the  Kolmogorov-Obukhov  theory  at  fre¬ 
quencies  aoove  1  GHz. 

Combining  the  general  expression  for  as  with  the 
transmission  equation,  we  get 


C  /// 
V 


GT(r)  GR(r) 
«£(£)  R^Cr) 


(3.18) 


where 


C 


3-m  ^2-m 


and  V  is  the  total  scattering  volume. 

In  the  SIGNATRON  computer  prediction  model,  the  path  loss 
(3.18)  is  calculated  using  a  numerical  integration  procedure. 

This  procedure  accurately  includes  the  effects  of  the  actual  beam 
pattern,  including  sidelobes.  The  computer  program  does  this 
simultaneously  for  several  receiver  beams  so  that  diversity  per¬ 
formance  can  be  estimated  accurately.  The  program  also  models 
long  term  variability  adapted  from  the  NBS  method.  Climatological 
effects  must  be  taken  into  account  when  selecting  the  atmospheric 
parameters,  and  LQ.  The  path  loss  computed  by  the  numerical 
integration  procedure  contains  contributions  from  both  the  basic 
path  loss  and  the  aperature-to-medium  coupling  loss.  This  result 
is  used  in  all  subsequent  calculations  for  digital  modem  pre¬ 
dictions.  In  addition,  the  closed  form  approximations  previously 
described  are  computed  as  a  means  of  verifying  the  utility  of 
these  expressions  for  applications  where  numerical  integration  is 
not  possible  or  convenient.  In  this  calculation  the  path  loss  is 
evaluated  in  two  steps.  First,  the  basic  path  loss  with  wide  beams 
is  evaluated.  Then  the  increase  in  the  path  loss  due  to  narrow 
beams  is  found.  This  loss  is  commonly  called  the  aperture-to- 
medium  coupling  loss.  The  basic  path  loss  is  found  to  be 


4m  *  *I#  l0«  [°rS  *3'“  k2"* 


where 


*0 


k 

■ 

0 

d 


l 

■  variance  of  refractive  index  fluctuations 
(5.10-H  in  dry  air 

•  correlation  distance  of  refractive  index 
(typically  70  a) 

«  2*/i  ■  wavenumber  of  incident  wave 
«  spectrum  slope  (m«ll/3) 

■  scattering  angle 

«  distance  between  terminals 


This  simple  expression  is  valid  for  a  wide  range  topographical 
condition.  The  aperture-to-medium  coupling  loss  can  be  evalu- 
a ted  approximately  by  the  formula" 


L,  «  -10  1C  Vl  (*?  >  «?) 

(X  ♦  fm+lHa-J)'  ) 
^1  +  i  |  2(hR~hT)  +  © 


+10  log  [1  + 


+10  log 


RG  "  9RG 


I) 


where 


eTVf  0RV 

PH 


vertical  beam  width  of  transmitter  and 
receiver ,  respectively 

horizontal  beam  width  of  transmitter  and 
receiver 


p  <x,y)  »  l-d+x)1"®  -  (l+y)1’®  +  (l+x+y)1'® 


This  result  assumes  terminal  antennas  with  the  same  horizontal 
dimensions.  Results  for  the  nonsymmetric  case  have  been  pre¬ 
viously  presented. 


/ 


* 
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(3)  Path  Loss  Dependence  on  Key  Parameters 
The  total  path  loss,  L^>  is 


lT  ■  hbS  +  Lc'  G 

where  G  is  the  total  antenna  gain.  For  two  terminals  with 
antenna  diameter  0,  we  have 


G  it:  15  dB  +  40  log  +  40  log  * 

If  the  frequency  is  doubled  and  the  antenna  diameters  halved, 
then  G  and  Lg  are  unchanged.  From  the  expression  for  t  ^ 
tnen  have  the  simple  rule  of  thumb 

e  The  path  loss  increases  5  dB  for  each  factor  of 
two  up  in  frequency  if  antenna  dimensions  are 
decreased  proportionally. 

Hence,  existing  5  GHz  links  with  8*  antennas  would  perform  5  dB 
worse  at  10  GHz  with  4'  antennas  and  10  dB  worse  at  20  GHz  with 
2*  antennas.  The  attenuation  due  to  water  vapor  absorption  and 
rain  increases  these  losses  substantially  at  higher  frequencies. 

The  dependence  on  distance  is: 


e  The  path  loss  increases  3  dB  for  each  factor  of 
two  in  distance  as  long  as  the  scattering  angle 
is  constant. 

e  For  zero  elevation  angles,  the  path  loss  in¬ 
creases  8  dB  for  each  distance  octave  when  the 
coupling  loss  is  negligible. 


The  dependence  on  the  scattering  angle  is: 


e  The  basic  path  loss  increases  5  dB  for  each 

factor  of  two  in  the  scattering  angle  8,  assuming 
constant  path  length. 


e  The  aperture-to-medium  coupling  increases  up  to 
an  asymptotic  limit  of  9  dB  for  each  factor  of 
two  in  scattering  angle. 
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SECTION  4 
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PERFORMANCE  MODEL  FOR  DIGITAL  TROPOSCATTER 


Performance  of  a  digital  troposcatter  system  can  be  ex¬ 
pressed  in  terms  of  average  bit  error  rate  and  outage  prob¬ 
ability.  Although  the  average  BER  can  be  a  useful  parameter 
and  is  frequently  measured  in  experimental  tests,  it  lacks  a 
direct  relation  of  digitized  voice  quality  in  a  troposcatter 
application.  The  outage  probability  is  a  measure  of  the 
fraction  of  time  a  fixed  BER  threshold  is  exceeded,  i.e.  the  fraction 
of  time  the  system  is  "out".  The  outage  probability  represents 
one  point  on  the  cumulative  BER  distribution  function.  Defense 

Communication  System  (DCS)  standards  for  digital  links  are 

-4 

expressed  in  terms  of  a  7.5  x  10  outage  probability  with  a  BER  threshold 
-4 

criterion  of  10  ,  since  this  level  corresponds  to  the  start  of  HograHoH 

digitized  voice  performance  in  a  PCM  system.  In  this  section  we  develop 
an  extensive  performance  model  of  a  digital  troposcatter  sys¬ 
tem  which  provides  the  means  for  computation  of  bit  error  rate 
statistics. 

a.  Bit  Error  Rate  Statistics 

The  bit  error  rate  performance  of  a  digital  data  system 
operating  over  an  angle  diversity  troposcatter  link  is  dependent 
on  the  received  power  at  each  feedhorn,  the  correlation  between 
diversity  channels,  and  the  multipath  characteristics  of  the 
>  composite  diversity  system.  In  this  analysis,  we  develop  the 

mathematical  model  to  predict  (1)  the  bit  error  rate  statistics 
as  a  function  of  the  hourly  median  received  power  for  the  main 
beam  diversity  channel  and  (2)  the  long  term  probability  that 
•  the  bit  error  rate  is  exceeded.  Statistics  for  the  first  pre¬ 

diction  ere  based  on  the  short  term  (less  than  a  minute)  fading 
characteristics  of  the  troposcatter  channel.  The  long  term 
calculation  deals  with  the  availability  of  the  digital  tropo- 


4-1 


scatter  system  as  a  result  of  hourly  and  daily  variations  in 
path  parameters. 

Computation  of  bit  error  rate  statistics  implies  a  choice 
of  digital  modulation  format  and  receiver  processing  technique. 

We  have  selected  QPSK  for  the  modulation  format  and  a  general¬ 
ized  adaptive  decision-feedback  equalizer  (DFE)  for  the  pro-  ^ 

cessor.  QPSK  is  the  present  choice  of  the  two  developed  tropo- 
scatter  modems,  MD-918  and  the  DAR-4,  because  of  its  bandwidth 
efficiency,  nearly  constant  envelope  characteristic,  and  optimal 
detection  capability.  The  DFE  has  been  shown  in  extensive  field 
testing  of  the  MD-918  modem  to  provide  excellent  performance 
over  an  extremely  wide  range  of  channel  conditions  such  as 
might  be  anticipated  in  an  angle  diversity  application.  The 
MD-918  DFE  modem  was  also  the  modem  used  in  the  AAC  field  tests. 

The  general  problem  of  computing  the  short  term  bit  error 
rate  (BER)  distribution  for  a  fading  multipath  channel  is  a  dif¬ 
ficult  if  not  impossible  task.  One  can  determine  the  BER  dis¬ 
tribution  for  a  flat  fading  channel  with  diversity  order  D  in  a 
straightforward  manner  [4.1]  but  the  introduction  of  implicit 
diversity  and  intersymbol  interference  effects  due  to  multipath 
precludes  a  closed  form  transformation  of  variables.  The 
short  term  average  BER  for  a  space  or  frequency  diversity  sys¬ 
tem  has  been  determined  following  the  procedure  developed  in 
[4.2]  with  a  resulting  expression  in  the  form  of  a  determinant. 

The  BER  distribution  can  be  computed  from  knowledge  of  the 
eigenvalue  structure  of  the  implicit  diversity  configuration 
that  results  from  this  procedure.  Two  important  assumptions 
are  introduced  which  make  calculation  of  this  eigenvalue  struc¬ 
ture  possible.  The  calculated  eigenvalues  represent  a  decom¬ 
position  of  the  troposcatter  channel  into  independent  implicit  * 

diversity  channels  with  average  received  powers  equal  to  the 
eigenvalue  set.  This  transformation  to  an  equivalent  independent 
diversity  system  then  allows  the  computation  of  the  mean  BER  and 
the  BER  distribution.  In  the  sequel,  the  analysis  developed  in  * 
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[4.2]  is  extended  to  include  angle  diversity  configurations 
and  calculation  of  the  BER  distribution. 

(1)  Communication  System  Definition 

The  communication  system  model  under  consideration  is  shown 
in  Fig.  4.1.  Complex  notation  is  used  to  represent  in-phase  and 
quadrature  components  and  explicit  modulation/demodulation  opera¬ 
tions  are  not  shown.  The  data  to  be  transmitted  (s^)  are  selected 
from  the  set  (+l+j)  for  a  quaternary  phase-shift-keying  system. 

The  transmitted  waveform  is 

oo 

s(t)  =k=£00  Skf  (t-kT)  (4.1) 


where  T  is  the  interval  between  successive  symbol  transmission, 
i.e.,  symbol  interval.  The  data  rate  for  this  QPSK  system  is 
2/T  bits/second.  The  transmitter  impulse  response  f(t)  is  de¬ 
fined  as  a  unit  energy  real  function. 

The  fading  multipath  channel  can  be  represented  by  an 
ensemble  of  zero-mean  random  functions  with  complex  Gaussian 
statistics.  Individual  diversity  channels  in  an  angle  diver¬ 
sity  system  will  exhibit  correlation  between  each  other  which 
can  be  completely  described  by  the  second  moment  of  the  chan¬ 
nel  impulse  response  h^(t), 

hi(t)  h^(x)  =  6(t-x)  Efa  j  “  1,2, . .  ,D, 

(4.2) 

where  6(t)  is  the  impulse  function,  E^  is  the  mean  received 
energy  per  bit  for  the  main  beam  diversity  channel,  rep¬ 
resents  a  relative  received  power  difference  between  the  diver¬ 
sity  channels,  and  p^j(t)  is  a  unit  area  function  equal  to  the 
mean  squared  envelope  of  the  impulse  response.  It  represents 
the  average  cross  power  response  at  each  delay  value  to  inrpulse 
function  excitation  and  will  be  referred  to  as  the  cross  multi- 


Fig.  4.1  Communication  System  Model 


path  profile.  Twice  the  RMS  width  of  the  profile  is  the  cross 
multipath  spread#  2c^j.  For  convenience  the  main  beam  diver¬ 
sity  channel  is  designated  as  channel  1  which  requires  A^  =  1. 
In  a  typical  angle  diversity  application  utilizing  dual  vertical 
feedhorns  on  each  antenna  with  odd  numbers  designating  main 
beam  diversities  and  even  numbers  designating  elevated  beam 
diversities,  the  characterization  has  the  special  form 


Aii  =  1  '  pii^  *  1  odd  (4.3.1) 

pii(t)  =  p2(t);  i  even  (4.3.2) 

pi,i+l  *  p12 ^ '  i  odd  (4.3.3) 

pij(t)  *  0  •  1  ^  3  and  i  (odd)  fi  j-1  (4.3.4) 


Conditions  (4.3.1)  and  (4.3.4)  state  that  all  main  beam 
diversities  have  equal  statistics  and  are  independent  of  each 
other.  Conditions  (4.3.2)  and  (4.3.4)  state  that  all  elevated 
beams  have  equal  statistics  and  are  independent  of  each  other. 
Finally,  conditions  (4.3.3)  and  (4.3.4)  state  that  the  cross 
multipath  profile  is  the  same  for  each  feedhorn  pair  but  that 
the  correlation  between  feedhorn  ports  on  different  antennas  is 
zero.  The  elevated  beam  powers  A..  .  even,  in  a  vertical  angle 

li  t  X 

diversity  system  generally  have  less  received  power  because  of 
the  larger  scattering  angle,  but  nonhomogeneous  conditions  can 
cause  A^  to  be  greater  than  unity  some  fraction  of  the  time. 

The  definition  of  as  the  received  energy  per  bit  for  the 
main  beam  diversity  A^1  =  1,  follows  from  the  fact  that  one  re¬ 
ceived  bit  is  represented  by  the  pulse 


00 

+  1  /  f(t-u)  h.(u)  du 

—00 


which  has  energy 


J  dtl  J  f  (t-u)h^(u)du|*  -  Au  J  dt  J  f^t-u)  ^P^u)  <*« 

m  \  I.Pl(tt)  dn  *  K 

The  ensemble  representation  is  used  here  because  the 
analysis  to  follow  assumes  a  slowly  varying  multipath  channel 
which  the  equalizer  can  track.  Hence  the  BER  statistical  per¬ 
formance  is  determined  by  computing  the  channel  ensemble  sta¬ 
tistics  of  the  bit  error  probability  for  the  Decision-Feedback 
Equalizer  (DFE) .  An  overbar  is  used  to  denote  channel  ensemble 
averages.  Brackets  <>  will  be  used  to  denote  average  over 
noise  or  source  statistics  for  a  particular  channel  realization. 
For  the  white  Gaussian  noise  channels  the  additive  noise  terms 
n^(t)  are  zero  mean  and  have  second  moments 

(Bj (t)  Bj  (t))  * 

where  Nq  is  the  noise  spectral  density  in  watts/Hz. 

Each  forward  filter  in  a  realizable  DFE  consists  of  a 
finite  length  tapped  delay  line  filter  with  tap  spacing  of  t 

s 

seconds.  If  there  are  K2  "late"  taps  and  "early  taps",  the 
tapped  delay  line  filter  impulse  response  is 

**  . 

w,  (t)  -  £  v  6  (t-kr  )  (4.4) 

1  k- ^  **  * 

where  the  tap  gain  value  is  chosen  with  a  complex  conjugate  for 
notational  convenience  later.  The  receiver  filter  has  impulse 
response  f(-t)  which  matches  the  transmit  pulse  waveform.  For 
this  choice,  in  the  absence  of  multipath,  and  after  time  syn¬ 
chronization,  the  optimum  tap  weights  for  the  main  beam  diver¬ 
sity  forward  filtei:  reduce  to  one  tap  on  and  the  rest  off.  The 
output  of  the  forward  filter  is  sampled  at  the  symbol  rate  1/T 
to  produce  the  sample  xm.  The  backward  filter  correction  sample 
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is  summed  with  the  sampled  receiver  filter  output  x  to  pro- 

111 

duce  the  detection  variable  2  at  the  mth  symbol  spoch.  The 

m  /s 

backward  filter  weights  previous  decision  i>o,  by  the 

complex  weight  b^  to  form  the  sample  output 


B 

I 

i=l 


b.s  . 
1  m-1 


(4.5) 


at  the  decision  time  for  the  s  symbol.  The  decision  on  s  is 

m  m 

denoted  as  sffl.  The  decision  process  takes  the  form 


z  *  x  +  y 

m  m  m 

s  »  sgn  (Re (z  ) )  +  jsgn  (Im(z  )) 
m  m  in 

where  sgn  is  the  signum  function. 

The  parameters  of  the  DFE  are  the  number  of  forward  filter 
taps  K,  the  forward  filter  tap  spacing*  t#  and  the  number  of  back¬ 
ward  filter  taps*  B.  The  optimum  DFE  requires  IC  »  •  »  B  and 

■  1/to  where  »  is  the  two-sided  bandwidth.  A  practical  choice 
of  parameters  for  tropo scatter  channels  was  determined  from  com¬ 
puter  simulation  of  a  fading  channel  equalizer  application  to  be 
K  -  3  «  B  and  t#  -  T/2. 

' (2)  DFE  Performance  Analysis 

Previous  methods  of  determining  bit  error  rate  (BER)  per¬ 
formance  of  practical  equalizer  structures  have  been  restricted 
to  Monte  Carlo  simulations  £4.3,  4.4]  using  an  ensemble  of 
multipath  channels.  An  analytic  approach  for  calculating  the 
BER  is  complicated  by  the  presence  of  an  intersynibol  interference 
term  in  the  signal-to-noise  ratio  (SNR)  expression  for  a  partic¬ 
ular  channel  realization.  Because  of  this  term,  the  calcula— 


tion  of  the  average  over  the  channel  ensemble  is  a  formidable 
task.  On  the  other  hand,  omission  of  the  intersymbol  inter¬ 
ference  term  leads  to  a  mean  BER  which  for  a  fixed  SNR  mono- 
tonically  decreases  with  increasing  multipath  spread.  For  any 
practical  equalizer  under  a  fixed  SNR  constraint  the  mean  BER 
will  initially  decrease  for  increasing  multipath  spread  but 
then  increase  as  the  multipath  exceeds  the  equalizer's  capa¬ 
bility  to  mitigate  the  intersymbol  interference.  Elimination 
of  the  intersymbol  term  provides  a  convenient  lower  bound 
which  both  shows  the  intersymbol  interference  penalty  and  is 
an  accurate  performance  estimate  for  small  multipath  spreads. 
However,  the  lower  bound  is  too  loose  for  performance  calcu¬ 
lation  when  the  rms  multipath  spread  is  on  the  order  of  the 
forward  filter  width.  The  analytic  procedure  to  be  developed 
allows  accurate  calculation  of  the  lower  bound  for  no  inter¬ 
symbol  interference  and  by  an  approximation  of  the  intersymbol 
interference  effect  provides  a  performance  estimate  for  the 
large  multipath  case. 

Two  key  assumptions  are  introduced  in  the  analysis  of  the 
performance  of  a  DFE  on  a  fading  complex  Gaussian  channel.  These 
assumptions  have  been  used  previously  [4.2,  4.5]  to  derive  predic¬ 
tions  of  average  bit  error  rate  (ABER)  for  conventional,  i.e., 
2S,  2F,  and  2S/2F,  diversity  configurations.  Test  results  both 
on  a  channel  simulator  and  in  field  test  measurements  showed 
excellent  agreement.  These  assumptions  are  again  used  and  the 
analysis  extended  to  include  the  prediction  of  ABER  for  angle 
diversity  configurations  and  the  prediction  of  the  BER  distri¬ 
bution  for  both  conventional  and  angle  diversity  systems.  The 
BER  distribution  predictions  will  be  used  subsequently  to  de¬ 
velop  an  outage  probability  prediction  on  a  yearly  basis. 

We  begin  the  analysis  by  considering  a  general  Xth  order 
diversity  system.  In  the  notation  of  the  previous  subsection, 
the  received  signal  for  the  ith  diversity  path  has  the  form 


fA(t)  ■  £  «a /h ^  (ac)  f(t-mT-x)  dx  +  nA(t) 


i  ■  1,2, ...I 


(4.6) 


We  define  the  combined  transmitter  and  receiver  response  as 


g(t)  •  J  f(t+x)  f(x)  du 


and  its  convolution  with  the  channel  as 


SB 

•yflj’q^t)  *  J  g(t-u)  h^u)  du. 


(4.7) 


For  a  sampling  time*  t^r  i  *  1,2, ...I,  the  detection  variable  z 
has  the  form 

I  K. 
z 


£1  k~K  ["»  ■£-  *»  1 

where  the  K-Kj+Kj+l  forward  filter  taps  per  diversity  are 
apportioned  between  Kj  "past*  taps  and  Kx  "future"  taps  and  one 
center  tap.  The  noise  variate 

vi<ti)  -  yn^tj-u)  f(u)  du 

is  zero-mean  Gaussian  with  second  moment 


<v^(tj-JcT)  Vj  (tj-lT)>  ■  0  i^j  (4.8a) 

^(t^kT)  vj  (tj-lT)>  -  N0  g(kT-iT).  i-j  (4.8b) 


•because  the  mean  delay  through  individual  angle  diversity 
channels  is  different,  different  sampling  times  must  be 
assumed  for  generality. 
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Thera  exists  an  obvious  representation  in  a  K^Kj+Kj+l  dimen¬ 
sional  column  vector  form,  viz., 

2  -  £  (v£  JL  *-9-i  +  2i 

where  the  accent  mark  refers  to  complex  conjugate  transpose. 
The  vector  w^  represents  the  forward  filter  complex  tap  gains 
for  diversity  i,  represents  the  ith  diversity  sampled  con¬ 
tinuous  filter  response  for  the  mth  transmitted  symbol,  and  v^ 
is  a  zero-mean  Gaussian  noise  vector  process  with  positive 
definite  covariance  matrix  NqGq  where 


) 


(4.9) 


Gq (k, A)  -  g(kr-ix) 


(4.10) 


Por  QPSK,  every  other  bit  decision  is  made  on  the  real  part 
of  z  and  there  is  quadrature  symmetry  and  thus  the  real  noise 
power  affecting  that  decision  is 


<*e(£  w'v,)2>  -  X,  wiG  w. 
1=1  11  2  i-1  101 


(4.11) 


If  one  assumes  s  ■  0,  m  0,  there  is  no  intersymbol  inter- 

m 

ference  and  the  analysis  for  the  BER  statistics  leads  to  an  SNR 
expression  for  a  particular,  channel  realization  which  is  a  qua¬ 
dratic  form  of  the  type  g£  go1  3o  where  3q  is  a  conP3-®*  Gaussian 
vector  with  statistics  determined  by  the  channel  and  GQ  is  a 
positive  definite  matrix.  After  a  diagonalization  procedure  the 
bit  error  rate  statistics  as  a  function  of  the  SNR  can  be  com¬ 
puted  for  differentially  coherent  detection.  This  calculation 
leads  to  lower  bound  expressions.  When  the  intersymbol  inter¬ 
ference  is  present,  an  effective  signal-to-noise  ratio  can  be 
defined  which  leads  to  a  quadratic  form,  but  both  the  vector  and 
the  matrix  have  random  components.  A  method  of  deriving  sta- 
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tistics  of  this  quadratic  form  is  not  apparent.  One  course  which 
has  proved  effective  is  to  approximate  the  intersymbol  inter¬ 
ference  effect  in  a  manner  which  will  modify  the  matrix  in  the 
quadratic  form  in  order  to  keep  its  deterministic  nature.  This 
can  be  accomplished  if  the  following  two  approximations  are 
made. 

(1)  Assume  the  interference  symbols  a.  m  j*  0,  are 
Gaussian  distributed  rather  than  complex  binary. 

(2)  Approximate  the  power  from  the  fading  in  ter  symbol 
interference  by  its  mean. 

This  approach  assumes  that  the  intersymbol  interference 
(ISI)  after  equalization  in  a  fading  multipath  channel  can  be 

well  approximated  by  an  equivalent  non-fading  additive  Gaussian 
noise  term.  Since  after  equalization,  the  ISI  is  generally 
small,  most  bit  errors  are  a  result  of  fades  in  the  desired 
signal  component.  Thus,  approximation  of  the  interference  char¬ 
acteristics  does  not  significantly  alter  the  BER  statistical 
analysis. 

In  the  absence  of  error  propagation  which  can  be  neglected 
in  the  analysis,  the  backward  filter  cancels  the  interference 
contribution  due  to  past  symbols ,  i . e . ,  symbols  which  arrive 
before  the  symbol  currently  being  decided  on.  If  we  arbitrarily 
take  the  symbol  selected  for  decision  as  sQ,  the  remaining  inter¬ 
ference  is  classified  as  future  and  the  contribution  for  the  ith 
diversity  is  given  by 
M 

E  Vtfami  *  -i-i* 


where 

M 

-1  ” 

is  defined  as  the  ith  diversity  ISI  vector. 


(4.13) 
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The  sum  can  be  truncated  after  a  few  terms  as  the  dot  pro¬ 
duct  will  disappear  for  large  m  due  to  the  finite  duration  of 
the  pulse  function  g^ (t) .  The  first  approximation  is  to  take 

s  as  a  zero  mean  complex  Gaussian  variate  with  mean  magnitude 
m  2 

squared  value  \»  •  This  quantity  should  be  less  than  two  because 
the  approximation  of  a  unit  magnitude  binary  variable  by  a  unit 
variance  Gaussian  variable  will  certainly  lead  to  pessimistic 
performance  results  due  to  the  tails  of  the  Gaussian  distribution. 
For  large  multipath  spreads,  a  choice  of  v  *1/2  was  found  to 

provide  good  agreement  between  calculated  and  measured  values 
in  a  average  BER  evaluation  [4.2,  4.5]  for  conventional  diversity 
systems. 

Equation  (4.9)  can  be  rewritten  as 
I 

2  ■  £.  (aoi^b  S°  +  -i  +  Hi)  <4-14> 


where  the  desired  signal  component,  the  noise  component,  and  the 
ISI  component  are  shown  explicitly.  Note  that  the  power  in  the 
signal  component  fades  due  to  random  variation  in  the  vector 
and  that  the  power  in  the  ISI  component  also  fades  due  to 
the  dependence  on  the  random  vectors  q^ ,  m  **  1,2,...M.  The 
noise  power  component  does  not  fade  as  has  been  shown  in  Eq. 
(4.11)  where  a  fixed  noise  power  independent  of  the  transmission 
medium  variations  results.  The  second  assumption  leading  to  a 
tractable  analysis  is  to  assume  that  the  power  from  the  fading 
ISI  components  can  be  approximated  by  its  mean.  Let  6^  repre¬ 
sent  the  actual  detected  power  from  the  noise  and  ISI  components. 
Since  the  additive  noise  contribution  from  v^  is  independent  of 
the  source  digit  contribution  from  u^,  V  i,  j,  the  effective  noise 
power  for  the  decision  on  the  real  part  of  z  is 

«2  -  <[r*(£i  wi(vi  +  Ui,)]> 


(4.15) 


Note  that  in  forming  the  real  part  in  (4.15)  for  complex 
representation  of  bandpass  waveforms#  terms  of  the  form 

-i-i-j-j 

are  double  frequency  terms  and  average  to  zero.  The  1SI  co- 
variance  matrix  in  (4.15)  after  averaging  over  the  source 
statistics  is 
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We  now  find  it  notationally  convenient  to  define  the  ex¬ 
tended  ik  column  vectors 


r-1! 

r  9mi  I 

• 

• 

e 

'  S*  - 

3m2 

e 

e 

• 

.Wj  . 

-  3ml  - 

the  K  x  K  sub-matrices 


and  the  IK  x  IK  extended  matrix 


Go  + 

S11 

S12  *** 

S1X 

S21 

G0  + 

S22  S23 

...  S2I 

Sn  •••  ®X,X-1  G0  +  SII 


(4.19a) 


(4.19b) 


(4.19c) 


The  matrix  GQ  has  been  previously  defined  in  (4.10).  The 
outer  product  c^.g^  is  related  to  the  cross  channel  multipath 
profile  channel  moments 


h^tjhjtt)  «  «(t  -  T>*b*J.jPj.j*t)  *  *  1,2, ...I 

(4.20) 
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We  can  write  the  outer  product  contributions  to  (4.19b)  as 
a  K  x  K  matrix  with  elemental  values  at  the  k,£  location  equal  to 


*  qi(ti-mT-kTs)qJ(t^-xnT-4T#)  ,  i,j  -  1,2,... I 

(4.21a) 

Performing  the  indicated  integration  in  Eq.  (4.7)  and  using 
the  second  moment  result  of  Eq.  (4.20),  we  have 

(m)  »  A^  J g(t^-«(T“hTa-x)g(tj-nff“iT#-x)pij  (x)dx 

•  « 

(4.21b) 


This  set  of  integrals  (Eq.  4.21b)  completely  describes  the 
transmission  medium  effects  on  the  DFE  performance. 

The  noise  and  ISZ  power  in  the  real  channel  symbol  detector 
reduces  from  (4.18)  to 


S2  . 

2  ~  ~ 


(4.22a) 


and  the  signal  power  at  the  same  decision  device  is 


[nelw'aol] 


(4.22b) 


2 

We  define  signal-to-noise  ratio  (SNR)  as  Pg/2&  at  the 
detector  as  this  definition  is  consistent  with  the  SNR  definition 
which  provides  a  BER  of 


9ER1DEM,  ■  2'  *rf°  <VSN*)  (4.23) 

for  the  ideal  non-fading,  Gaussian-noise,  binary-symbol  channel. 
The  signal-to-noise  ratio  for  this  fading  channel  application  is 
then 
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p,5*  r»<s,’o>2/<!£,«>  14.24) 

which  we  seek  to  maximize  as  a  function  of  the  forward  filter  4 

weight  vector,  w.  Consider  the  generalized  dot  product  defined 
on  the  positive  definite  matrix  Gt  i.e. , 


(u,  v)  *  u'Gv. 

Equation  (4.24)  can  be  written  as 

a  [if  (w.c-1^)]2 

K0  (w#w) 


(4.25) 


(4.26) 


and  by  a  generalization  of  the  Schwartz  Inequality  we  obtain 

Eb  U«.G-1a0)|*  i  i  C4.271 

Pii£ — “•*  30 


with  equality  if  and  only  if 


^pt 


(4.28) 


Note  (4.28)  results  in  a  pure  real  term  in  (4.26)  to  satisfy 
the  Be  (*)  requirement. 

The  maximum  SNR  is  then 


.  ,  fb  r-l_ 

3max  „  *0  9o 

N0 


(4.29) 


* 
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As  a  check  on  the  normalization,  if  X  *  1,  a^«0  and  t^*0 
we  have  from  Eq.  (  4.7) 


q(-kT)  -  gt-k-r)  »  G- 


(4.3 


and  then 


~  g  (0)  -  ^ 
N0  N0 


(4.3 


as  is  expected  in  the  absence  of  . multipath  and  diversity.  In 

.  •  . .  . 

this  trivial  example  the  ideal  BEK  (4.23)  is  achieved  for  every 
fading  channel  realization. 

He  will  now  establish  that  the  mean  square  error  function 


<e2>  -  f  <1*  -  *>0|2> 


(4.3; 


which  is  minimized  by  the  decision-feedback  equalizer  leads  to 
the  same  set  of  forward  filter  weights  as  maximization  of  the 
SNR  in  (4.26)  for  the  appropriate  choice  of  scale  factor  d. 
After  squaring  and  averaging,  equation  (4.32)  becomes 

<*2>  *  ^'Soajx  -  ^  d  **(x'ao>  +  |  Ho£,Q£  +  *2 


-  I  \  !£'So-*r  +  No— ,GH- 


(4.3: 


Minimization  of  (4.33)  is  equivalent  to  the  LaGrange  multiplier 
problem 


I 


which  has  solution 


w  «  SG"1^.  (4.35) 

Since  Bq.  (4.28)  and  (4.35)  are  the  sane  except  for  an  unimpor¬ 
tant  scale  factor,  we  have  established  that  the  maximum  signal- 
to-noise  ratio  for  a  minimum  mean-square  error  DFE  and  a  particu¬ 
lar  channel  realization  is  the.  quadratic  form  given  by  Bq.  (4.29). 

The  maximum  SMB  result  (Bq.  4.29)  indicates  the  optimum  DFE 
performance  for  a  particular  fading  channel  realization.  Since 
is  a  function  of  the  transmission  medium,  the  BER  statistics 
must  treat  ^  as  a  random  process.  Note  our  two  approximations 
have  led  to  a  deterministic  form  for  the  matrix  6  which  makes  an 
analytic  treatment  feasible.  We  now  turn  our  attention  to  the 
evaluation  of  BER  statistics;  in  particular,  the  average  BER  and 
the  BER  probability  density  function.  We  begin  with  the  easier 
task. 

(3)  Average  Bit  Error  Rate 

The  average  bit  error  rate  is  computed  by  taking  the  average 
over  the  channel  ensemble.  The  bit  error  probability  for  QPSX 
for  any  member  of  the  channel  ensemble  is 

Pc  -  j  P*U  >  1)  +  ~  pr  (C  <  -1}  (4.36) 

where  C  is  a  zero-mean  Gaussian  random  variable  with  standard 
deviation  equal  to  /l/2pmJ1~.  The  bit  error  probability  is 


Pc  »  exp(-u2/2)du  *  i  erfc  /p—  (4.37) 

2p 

max 

The  subscript  c  denotes  that  the  detection  process  was  co¬ 
herent,  i.e.  the  receiver  knows  the  transmitted  carrier  phase. 
Differential  detection  of  phase-shift-keyed  signals  yields  a  bit 
error  probability  [4.6] 


Pd 


1 

2 


-P. 


max 


(4.38) 


The  tap  gains  of  an  adaptive  MMSE  equalizer  remove  the  phase 
and  frequency  difference  between  the  transmitter  and  receiver 
carrier  clocks  toan  ambiguity  of  180°  for  binary  transmission 
and  90°  for  quaternary  transmission.  This  ambiguity  is  eliminated 
by  differentially  encoding  the  transmitted  data.  The  performance 
of  coherently  detected  and  differentially  encoded  data  signals 
for  a  fixed  is  approximately  double  the  error  rate  given  by 

(4.37}.  When  practical  degradations  are  included  in  the  analysis 
the  performance  is  closer  to  the  differential  detection  expres¬ 
sion  (4.38).  Thus  the  ensemble  average  of  pQ  is  of  interest  as 
a  fading  channel  performance  bound  and  the  ensemble  average  of 
Pg  is  of  interest  as  an  estimate  of  realizable  modem  performance. 

Let  y(x)  be  the  probability  density  function  for  pm|)T  and 
Y(S)  its  Laplace  transform.  We  want  to  calculate  the  ensemble 
averages 


erfc  /x  y(x)dx 


(.4.39) 


0 


(4.40) 


The  Laplace  transform  Y(s )'  is  easily  obtained  after  a 
diagonal!  sat  ion  of  the  quadratic  form  in  Sq.  (4.29) 
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4 * .  „  4.. 

0  o 


(4.41) 


where  a  is  a  zero-mean  complex  Gaussian  vector  with  a  diagonal 
covariance  matrix  rij  *  i*l»2,...Ik,  where  I  is  the 

number  of  explicit  diversity  channels  and  k  is  the  number  of 
forward  filter  weights.  The  diagonalization  resulted  from  the 


transformation 


«  -  MG'^ 


(4.42) 


and  M  is  the  normalized  orthogonal  matrix  for  the  matrix 


G~**  gjjgj  G-1*,  i.e. 

G'*  3jlo  G‘%M  ■  MF 


(4.43) 


(4.44) 


M'M  -  X. 


(4.45) 


Thus  the  eigenvalues  X^  are  also  the  eigenvalues  of  the 

symmetric  matrix  G”*1  gflg^  G-**  and  the  unsymmetric  matrix 

G-*  g^jgj.  Since  the  components  of  a  are  uncorrelated  Gaussian, 
the  probability  density  for  Pinax  i»  the  convolution  of  exponen¬ 
tial  densities  which  yields  a  Laplace  transform  in  product  fora. 


y<8>  -  n  a  +  x.s) 

k-1  No  * 


(4.46) 


mean  BER  for  DPSX  is 


?d  -  |  Y(l)  -  i  XI  V 

2  2  k-1  * 


(4.47) 
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For  coherent  detection  it  is  necessary  to  find  the  coeffi¬ 
cients  in  a  partial  fraction  expansion  of  Y(s),  i.e.. 


Y(S) 


A  recursive  method  for  finding  the  coefficients  and  the  resulting 
error  probability  for  coherent  detection  is  detailed  in  [4.7] 
where  the  mean  BEK  for  the  infinite  length  optimum  single  pulse 
receiver  was. determined.  We  will  concentrate  on  the  DPSX  ex¬ 
pression  (4.47)  because  the  DPSX  bit  error  probability  curve  for 
the  Gaussian  noise  non  fading  channel  closely  matches  the  BEK 
for  the  HD-918  modem*  i.e. *  the  implementation  loss  and  differen¬ 
tial  encoding  loss  in  the  CPSK  detection  system  is  approximately 
equal  to  the  loss  incurred  in  ideal  DPSX  relative  to  ideal  CPSK. 
Figure  4.2  illustrates  this  result.  Since  in  a  fading  channel 
application  most  of  the  bit  errors  will  occur  during  periods  of 
low  SNR,  good  correspondence  between  ideal  DPSK  and  the  measured 
DFE  modem  at  low  SNR  is  desirable  for  a  DPSK  representation. 

For  a  DPSX  system  or  for  the  use  of  the  DPSK  expression  to 
approximate  the  Gaussian  noise  non  fading  performance  of  a  practi¬ 
cal  modem*  calculation  of  the  average  bit  error  probability  does 
not  require  the  determination  of  the  eigenvalues.  Mote  that 
(4.47)  has  the  determinant  form 

Pa  -  i  M.t|I  +  G^Cl)'1  (4.48) 


where  the  covariance  matrix  C  is  defined  by 


(4.49a) 


BIT  ERROR  PROBABILITY 


The  matrix  C  is  of  rank  IK  and  has  a  submatrix  structure 
of  the  form 


R11  R12 


*21  *22 


*11  *I-lf  2 


•  •  •  Ri 


(4.49b) 


where  R^  is  a  K  x  K  matrix  defined  by 


*ij  *  *0iq0j 


(4.49c) 


and  thus  it  has  elemental  values  at  the  k, t  location  equal  to 
*£(0)  which  was  defined  in  Eq.  4.21b. 

The  general  result  (Eq.  4.48)  is  valid  for  any  Ith  order 
diversity  system.  In  subsection  4. 1.2. 3  we  examine  how  these 
matrices  simplify  when  specific  diversity  configurations  are 
assumed.  First  we  derive  the  formulas  for  the  BER  probability 
density  function  and  describe  how  a  yearly  outage  probability 
estimate  can  be  obtained. 


Bit  Error  Rate  PDF  and  Outaae  Probabili 


For  a  particular  channel  realization  the  signal-to-noise 
ratio  has  been  expressed  as  the  sum  (c.f.  Eq.  4.41), 


■  i  i  • 


(4.50 


The  random  fading  causes  the  to  be  independent  zero-mean 
complex  Gaussian  values  with  variance 


i  »  1,2. ..IK 


(4.51) 


where  the  X^  are  eigenvalues  of  the  non  symmetric  matrix  G~XC. 
The  matrices  G  and  C  are  defined  by  Eqs.  (4.19)  and  (4.49), 
respectively.  The  bit-error-rate  is  itself  a  randan  process 
equal  to 


i  ~P 
1  _  max 


(4.52) 


The  probability  density  function  (PDF)  of  pfl  can  be  found 
by  a  transformation  from  tne  PDF  of  P _ .  The  random  variable 

i  i2 

|a^|  has  an  exponential  PDF  whose  Laplace  transform  is 


Yi(s)  -  (1  +  X±s) 


(4.53) 


and  since  theja^p  are  independent,  the  Laplace  transform  of  the 


PDF  for  is 


-  IK  E.  . 

Y(s)  -  Jy(p)e~spdp  -  H  (1  +  -2  XjS)*1 

0  i-1  N0 


(4.54) 


The  X^  values  are  the  eigenvalue  set  from  the  matrix  g“aC, 

Partial  fraction  expansion  techniques  discussed  in  the 
previous  subsection  can  be  used  in  a  straightforward  but  some¬ 
what  tedious  calculation  to  obtain  the  PDF  y(p).  To  calculate 
the  BER  PDF  x(p)  we  use  (4.52)  in  the  transformation 


x(p)  -  y(p)  |j£|l  ■  i  y  (ln(l/2p) ) 


(4.55) 


p  -  tn(l/2p) 


The  outage  probability  has  been  defined  as  the  probability 
that  the  bit  error  rate  exceeds  a  value  p  and  is  equal  to  the 
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cumulative  distribution  of  the  PDF.  This  distribution  is  a 
function  of  the  mean  bit  energy  to  noise  power  ratio 

V  -  VHo 

or  in  dB  units  r  *  10  Log(y)  dB. 

Thus ,  we  have  the  conditional  outage  probability 

B(p|D  *  prob{BER  £  y)  * 

P 

The  integration  in  (3.56)  is  more  easily  accomplished  before 
the  transformation#  i.e. 

B(p|D  -  PrGb(PinaX  <  4n(l/2p)|^/ta0  -  y) 

4n  (l/2p) 

y(p)<ap.  (4.57)’ 

0 


The  unconditional  outage  probability  averages  over  the 
long  term  variations  in  which  occur  on  any  fading  channel. 

In  troposcatter  propagation  these  variations  over  a  year  are 
well  described  by  the  Gaussian  PDF 


f(T) 


1  e-QT-m)2/2o2 

a/2  rr 


•  <!*<• 


(4.58) 


where  m  and  o  are  the  mean  and  standard  deviation  of  the  long 
term  variations  expressed  in  dB  units.  Ate  unconditional 

outage  probability  for  troposcatter  is  then 


B(p) 


x(q|r)f(T)dr  dq 


kThis  calculation  assumes  that  the  multipath  distribution  can 
be  approximated  by  its  mean  value.  This  assumption  is  sug¬ 
gested  by  conditional  outage  probability  results  presented 
later. 
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where  the  dependence  of  the  BER  PDF  on  r  has  been  shown 
explicitly.  The  unconditional  outage  probability  can  be 
expressed  as  an  average  of  the  conditional  form  in  (4.56)  by 


B  (p) 


f(r)B(p|r)dr 


(4.59) 


.  In  practice  the  calculation  of  the  above  integral  is  diffi¬ 
cult  because  each  value  of  r  requires  a  different  partial  frac¬ 
tion  expansion  solution  to  obtain  x(p|r)«  Fortunately*  B(p|r) 
is  a  very  steep  function  of  r  relative  to  the  PDF  f(T)  so  a 
relatively  simple  numerical  integration  routine  should  approxi¬ 
mate  B (p)  for  a  few  fixed  values  of  p  quite  well.  Of  particular 

-2  -4  - 

interest  in  this  analysis  are  the  values  p  “  10  *  10  *  and  10 

A  numerical  procedure  to  estimate  B(p)  might  proceed  as  follows. 
First  find  the  value  of  SNR  from  the  analysis  of  the  previous 
subsection  which  gives  a  mean  BER  of  p.  Call  this  value  rQ. 

When  r  -  rQ  the  instantaneous  BER  is  most  of  the  time  much 
smaller  than  the  average  BER.  Hence*  non  negligible  values  of 
x(p|r)  must  be  for  a  T  value  smaller  than  Le t  the  SNR  . 

decrease  in  a  dB  steps*  i.e. 


ri+l  “  ri”&  1  "  O*1#2#*-* 


(4.60) 


and  find  the  corresponding  values  B(pjr^)  by  inverting  the 
relationship  given  in  (4.56).  This  involves  for  each  1*^  per¬ 
forming  the  partial  fraction  expansion  to  obtain  B(p|r)  and 
then  in  a  look-up  operation  determining  B(p|r)  for  a  specific 
p.  The  values  of  B(p|l*£)  should  start  small  and  increase  to¬ 
ward  unity.  After  Btpjr^)  exceeds  a  threshold  the  procedure 

can  be  terminated  and  B(p)  approximated  by 

N  • 

b(p)  4  S  f<ri)B(p|ri)  /  fflr^  (4.6D 

In  this  manner  the  outage  probability  distribution  for  the  year 
can  be  determined. 
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b.  Special  Diversity  Configurations 

In  this  subsection  we  consider  specific  diversity  configura 
tions  of  practical  interest  and  we  examine  what  simplifications 
in  the  previous  analysis  are  possible.  The  two  special  cases 
of  interest  are  identical/independent  diversity  systems  and 
angle  diversity  systems.  The  former  includes  conventional 
2S/2P  diversity  and  the  latter  includes  the  angle  diversity 
configurations  such  as  2S/2F/2A,  2S/2A,  and  2F/2A. 

(1)  Identical/Independent  Diversity 

For  this  example  if  the  number  of  identical/independent 
diversity  paths  is  D  then  the  eigenvalues  are  repeated  D  times. 
This  result  follows  from  the  identical  and  independent  statis¬ 
tical  nature  of  each  diversity  path*  i.e.  the  fading  channel 
impulse  response  is  characterized  by 


hA  (t)hj*  (t)  -  fttt-TjEj^^ptt),  i,j«l,2,...D 

The  set  of  integrals  (Bq.  4.21b)  which  describe  the  effects 
of  the  transmission  medium  on  the  DFE  performance  reduce  to 


^  (m) ^  g  (t^-«i*-i»i^-«»3 


.-mT-kT  -x)g  (t ,  -mT-4,T  -x)p  (x)dx 


*£[<m)  "  0  i#j 


For  statistically  identical  diversities  the  mean  path  delay 

through  each  of  the  diversity  paths  is  the  same.  The  sampling  times 

ii 

t^  are  all  chosen  equal  and  hence  (a)  does  not  depend  on 
!•  The  extension  to  fora  the  matrices  G  and  C?  (4.19c)  and 
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* 


(4.49b),  respectively,  thus  requires  the  submatrix  Gg+S^  in  G 
and  R^  in  C  to  be  repeated  on  the  diagonal  D  times.  Since 
the  eigenvalues  are  determined  from  the  matrix  G  *C,  the  K  dis¬ 
tinct  implicit  diversity  eigenvalues  are  computed  from  the 
matrix 


(VS11>  *11* 


For  this  special  case  we  can  summarize  the  procedure  for  eval¬ 
uating  performance  as  follows.  First  let 


R 


Rjj#  R(k#<t>) 


g  (tj-lcTg-xJg  (tj-i^-xjp  (x)dx 


s 


*11' 


S  Oc,  4) 


¥i 

”o  m*l 


) 


GQ(lc#t)  ■  g(lcta-4TB)« 

where  1c  and  t  represent  the  X  integer  values  indicating  the 
equalizer  tap  locations.  The  integer  M  is  the  number  of 
significant  future  symbol  interferers.  The  average  BER  is 
equal  to 

*1  '  2  («•*!*  +  ^  <G0+S) 

To  find  the  BER  PDF  and  the  outage  probability  the  pro¬ 
cedure  outlined  in  the  previous  subsection  is  to  be  followed 
using  the  positive  definite  eigenvalue  set  repeated 

D  times.  The  eigenvalue  set  is  determined  from  the  non  sym¬ 
metric  matrix  (0Q  +  S)”*R, 


9 
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» 


* 


X 


(2)  Angle  Diversity 

Consider  a  diversity  configuration  with  I  angle  diversity 
paths  and  D  independent  frequency  or  space  diversity  paths. 

The  effect  of  the  D  independent  diversities  is  to  repeat  the 
IK  eigenvalues  of  the  correlated  diversity  system.  Since  the 

extension  of  the  result  for  Et>l  is  clear,  we  restrict  ourselves 

« 

here  to  examining  the  pure  angle  diversity  system,  i.e.  0*1  • 

We  now  wish  to  establish  that  under  certain  conditions, 
the  correlation  between  an  angle  diversity  pair  does  not  signi¬ 
ficantly  effect  performance  and  it  can  be  neglected.  Under 
this  condition  angle  diversity  reduces  to  an  independent  diver¬ 
sity  example  but  with  different  diversity  branch  strengths  and 
multipath  profiles. 

Consider  dual  vertical  angle  diversity  (1*2)  with  a  K 
tap  forward  filter  in  the  DPE.  The  taps  are  equally  spaced 
with  a  width  of  T/2.  The  transmission  medium  is  completely 
characterized  by 

hj(t)Y(T)  -  Jlt-TjA^tt) 

h2(t)h2*(T)  ■  a  (t— r)A22P22  (t) 

hi<t)l»2*<T)  ■  a  <t-T)A12P12  (t) 

The  multipath  profiles  p^(t)  have  rms  multipath  spreads 
of  Oj,  where 

✓  ■  2 

a\mJ  t2pil(t)dt -Ij  t  pu(t)  dt  J  ,  i-1,2. 

0  '  0  ' 


4-29 


For  vertical  angle  diversity  Oj  is  normally  a  little 
larger  than  o^.  The  correlation  between  angle  diversity  chan¬ 
nels  for  CW  transmission  is 


SCW 


/*11*22 


|hl(t)dt  J h2* 


(s)ds 


‘12 


•  J  *11*22 


where  we  have  used  the  unit  area  definition  of  the  multipath 
profiles.  We  assume  that  %  is  not  negligible*  in  fact  for 
dual  vertical  angle  diversity  it  is  on  the  order  of  0.6.  cor¬ 
relation  effects  a  digital  transmission  system  through  the 
second  moment 


(m) 


g  (ti-mr-1CTs“X)g  (t^-mr-tTs-x)pij  (x) dx 


If  this  moment  is  negligible  when  i*j,  correlation  will 
not  effect  the  digital  system.  Note  that  the  autocorrelation 
function  g(*r)  has  width  -on  the  order  of.T  seconds*  .  NPP  be¬ 
cause  the  mean  path  delay  is  longer  for  the  elevated  beam*  we 
have  to  choose  a  later  sampling  time  for  the  elevated  beam 
diversity*  i.e.* 

H  -  *i  +  * 

where  &  is  the  estimate  of  additional  transmission  delay 
through  the  elevated  beam.  The  value  of  t  is  on  the  order  of 
2o^.  For  high  data  rate  systems  when  T  is  on  the  order  of  &* 
the  integral  d£J(m)  i«  small  compared  to  J  |  (m)  |  |^J(«)| 


x 


*  In  practice*  this  effect  is  schieved  by  adding  a  fixed 
delay  to  the  lower  beam  diversity. 
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except  when 

t,  +  A  -  lra  «w  t  -  k-r  # 

X  sis 

as  then  the  autocorrelation  functions  approximately  line  up. 

For  a  two  tap  forward  filter  (k*0,l)  with  ra“T/2,  the 
above  condition  can  not  be  satisfied  and  performance  is  approxi¬ 
mately  the  same  if  A^  wee  taken  to  be  zero.  When  an  additional 
tap  is  added  to  the  equalizer*  the  implicit  diversity  of  the 
system  increases  and  performance  improves.  The  improvement 
from  two  to  three  taps  is  not  particularly  large*  however. 

With  three  taps  (k*-l,Of  1)  there  will  be  some  off  axis  cor¬ 
relation  contributions  to  the  submatrices  when  k*l*  j*-l 
and  when  k*-l*  4*1.  Since  the  performance  improvement  between 
two  and  three  taps  is  not  large  and  since  the  two  tap  system  is 
not  affected  by  nonzero  correlation  >  0*  it  follows  that 
these  off-axis  submatrix  terms  can  be  neglected  for  a  three 
tap  DFE  evaluation.  Also  note  that  even  for  a  CW  system*  the 
'  neglect  of  a  correlation  equal  to  0.6  leads  to  only  about  a 
ldB  error. 

The  amission  of  the  off-axis  submatrix  terms  eases  the 
eigenvalue  determination  problem  because  with  independent 
angle  diversity  pairs*  the  three  tap  DFE  requires  the  calcula¬ 
tion  of  the  roots  of  two  cubic  equations  rather  than  calcula¬ 
tion  of  eigenvalues  from  a  matrix  of  rank  6. 

Zn  principle  these  considerations  suggest  that  the  cor¬ 
relation  effect  is  reduced  for  high  speed  data  systems  When 
the  data  symbol  interval  becomes  sufficiently  short  to  be  on 
the  order  of  the  mean  delay  between  angle  diversity  branches • 
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For  the  special  case  of  dual  vertical  angle  diversity 
with  neglected  correlation  between  the  angle  diversity  beams, 
the  performance  is  calculated  as  in  the  previous  subsection 
except  two  sets  of  K  eigenvalues  must  be  found  corresponding 
to  the  two  multipath  profiles  p^ (t)  and  p^ (t) .  These  eigen¬ 
values  are  weighted  by  the  angle  diversity  power  factors  A^ 
and . .  The  2K  eigenvalues  are  then  repeated  D  times  for  a 
D  identical/independent  configuration.  The  analysis  then  pro¬ 
ceeds  as  for  the  identical/independent  configuration  examined 
in  the  previous  subsection. 
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SECTION  5 


PERFORMANCE  STANDARDS  AND  EXAMPLE  COMPUTATION 


The  communication  channel  model,  path  propagation  prediction 
method,  and  modem  performance  calculation  technique  have  been  com¬ 
bined  for  prediction  of  digital  troposcatter  link  performance. 
Proposed  standards  for  U.S.  military  communication  networks  containing 
troposcatter  links  have  been  developed  by  the  Defense  Communica¬ 
tions  Engineering  Center  (DCEC) .  In  this  section  we  review  these 
standards  and  discuss  the  method  of  application.  The  RADC  tropo¬ 
scatter  test  link  is  then  used  as  an  example  for  comparison  of 
performance  prediction  against  the  DCEC  standard. 

a.  Performance  Standards 

Performance  standards  for  evaluation  of  U.S.  DCS 
digital  troposcatter  systems  are  based  on  the  communication 
requirement  established  in  DCEC  Technical  Report  12-76  dated 
November  1976.  This  report  identifies  five  ranges  of  outages 
and  performance  requirements  for  each  outage  range.  The  five 
outage  ranges  are  defined  as  follows: 

Range  I.  Outages  under  200  milliseconds  duration  which 
are  significantly  separated  in  time  from  any  other  outage. 
Such  outages  will  normally  occur  on  a  troposcatter  link 
operating  with  adequate  margin.  The  allowable  probability 
of  occurrence  of  these  outages  need  not  be  specified  since 
their  effect  on  the  user  is  trivial. 

Range  II.  Outages  with  durations  in  the  range  from  200 
milliseconds  to  five  seconds  which  are  significantly 
separated  in  time  from  any  other  outage.  Such  outages 
may  occur  on  either  a  line-of-sight  or  a  troposcatter 
link  operating  with  adequate  margin  (or  less  than  adequate 
margin  if  the  frequency  of  occurrence  is  high) .  The 
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allowable  probability  of  occurrence  of  such  an  outage 
should  be  specified  with  a  value  based  on  its  annoyance 
to  the  user. 


Range  III.  Outages  with  durations  greater  than  5  seconds 
but  less  than  2  minutes.  Such  outages  will  normally  occur 
on  a  line-of-sight  link.  The  allowable  probability  of 
occurrences  of  such  outages  should  be  specified  with  a 
requirement  more  stringent  than  that  for  Range  II  since 
such  an  outage  is  considered  to  be  a  disruption  to  the 
user  which  would  cause  call  termination. 

Range  IV.  A  recurrent  set  of  outages  of  average  duration 
under  200  milliseconds  which  occurs  at  an  average  rate 
of  from  2  to  5  per  minute.  This  condition  will  normally 
occur  on  a  troposcatter  channel  operating  with  barely 
adequate  margin.  The  allowable  probability  of  existence 
of  such  a  condition  should  be  specified  with  a  value 

-•M  - 

based  on  its  annoyance  to  the  user. 

Range  V.  Outages  with  duration  greater  than  2  minutes  or 
recurrent  outages  of  any  duration  which  occur  more  fre¬ 
quently  than  5  per  minute.  The  duration  of  any  such 
outages  or  period  of  recurrent  outages  should  be  included 
in  the  total  unavailability  specification  since  the 
channel  is  considered  essentially  unusable.  The  outages 
described  by  this  range  will  normally  occur  on  a  tropo¬ 
scatter  channel  with  long-term  fading  in  excess  of  the  fade 
margin,  or  on  a  line-of-sight  channel  experiencing  ano¬ 
malous  index  of  refraction  conditions. 
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The  troposcatter  channel  is  most  accurately  viewed  as 
having  a  compound  fading  characteristic  with  short  term  (e.g., 
hourly)  Rayleigh  fading  impressed  on  a  longer  diurnal  and  annual 
fading  distribution.  Historically,  the  performance  of  tropo¬ 
scatter  links  was  related  to  the  percentage  of  time  that  a 
particular  short  term  average  Signal-to-Noise  Ratio  (SNR)  could 
be  maintained  out  of  the  longer  term  fading  ensemble  (i.e., 
time  availability).  The  essence  of  this  convention  is  used 
in  DCEC  TR12-76  and  is  continued  here. 

The  short  term  average  SNR  is  defined  as  the  ratio  of 
the  average  received  bit  energy  divided  by  the  noise  power  in 
a  1  Hz  RF  bandwidth.  This  ratio  is  frequently  denoted  ^/Nq 
where  the  overbar  denotes  an  average  over  the  fading.  In  the 
comparison  of  diversity  configurations,  E^Ng  is  the  value 
associated  with  a  receiver  in  a  conventional  space/frequency 
diversity  system.  For  PSK  systems  E^  is  equal  to  the  average 
received  power  per  bit  times  the  symbol  duration,  i.e. 

*  Prt/L0G2  (M) 

where  P  *  average  received  power 

T  =  symbol  length 
M  *  number  of  bits  per  symbol. 

The  data  rate  in  bits/second  is 

R^  ■  LOG 2 (M) /T  bits/second 

and  thus  E^/Ng  can  be  expressed  as  the  ratio  of  average  re¬ 
ceived  power  to  noise  power  in  the  bit  rate  bandwidth,  i.e., 

V*o  -  VVb- 

Consistent  with  the  notation  in  TR12-76,  this  parameter,  short 
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term  SNR,  will  be  denoted  throughout  this  discussion  by  the 
symbol,  yq-  The  yq  necessary  to  obtain  a  specified  short  term 
performance  (here  described  on  the  basis  of  acceptable  outage 
rates  and  duration)  will  then  be  described  in  unavailability 
fashion  (i.e.,  to  be  maintained  for  all  but  a  certain  percent¬ 
age  of  time)  and  converted  into  required  system  gains  for  various 
DCS  digital  troposcatter  configurations.  The  performance  cri¬ 
teria  selected  in  TR12-76  to  be  used  in  these  determinations  are 
indicated  in  Table  5.1- 


Table  5.1 

DCS  PCM  VOICE  PERFORMANCE  CHARACTERISTICS 


A 

OUTAGE  RANGE 

_ 

B 

CRITERIA 

C 

OUTAGE  PROBABILITY  • 

II 

.2  sec  <:  outage  <  5  sec. 

7.5*  It)""4 

III 

5  sec  <  outage  <  1  min. 

7.5-10"5 

IV 

2  <  outages/min.  <  5 

outage  <  .2  sec. 

2.5* lO-3 

V 

5  <  outages/min. 

or 

outage  >  2  min. 

1-10"4 

An  analysis  in  TR12-76  relates  these  performance  criterion 
to  a  value  of  yq  exceeded  all  but  a  certain  fraction  of  time. 
Because  of  the  analytic  complexity  of  the  problem,  key  assump¬ 
tions  are  used  in  deriving  yQ  criteria.  These  assumptions  are: 

(1)  implicit  diversity  doubles  the  effective  order  of 
diversity  and 

(2)  the  outage  probability  as  a  function  of  Inr  is  taken 
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as  either  unity  or  zero  in  the  computation  of  yearly  outage 
probability. 


The  analysis  developed  here  extends  the  work  in  TR12-76 
by  calculating  the  exact  contribution  of  the  implicit  diversity 
and  the  exact  form  of  the  outage  probability  vs. SNR  for  inte¬ 
gration  over  the  yearly  SNR  distribution.  Significantly, 
these  new  results  show  that  implicit  diversity  does  not 
double  the  effective  diversity  order  and  that  the  piece-wise 
linear  approximation  of  the  outage  probability  is  poor  when 
comparing  diversity  configurations. 

Qualitatively  the  analysis  in  TR12-76  revealed  that  the 
L-band  (755-985MHz)  systems  were  moat  susceptible  to  outage 
range  condition  II  with  respect  to  the  performance  criteria 

in  Table  5.1.  This  suggests  an  L-band  yearly  outage  proba- 

-4 

bility  requirement  of  7.5  x  10  ,  i.e.,  the  probability  that 

the  HER  exceeds  10  4  for  0.2  sec.<  outage  <  5  sec.  over  a  year  is  less  than 
-4  ~ 

7.5  x  10  .  The  higher  frequency  C-band  systems  were  analyzed  in  TR12-76 

using  a,  5Hz  rms  Doppler  spread  in  the  evaluation.  As  a 
result,  the  outage  range  II  requirement  was  not  the  most  demand¬ 
ing  criterion  for  C-band  systems.  However,  measured  rms 
Doppler  spread  does  not  indicate  that  5Hz  spreads  will  be 
achieved  very  often  on  C-band  systems.  As  an  example. 

Figure  5.1  compares  the  measured  2a  rms  Doppler  spread 
for  a  168  mile  path  using  both  C  band  and  L  band  transmis¬ 
sions.  The  small  difference  between  these  distributions 
suggests  that  the  outage  range  II  criterion  should  be  applied 
to  all  troposcatter  systems  regardless  of  frequency. 

This  result  can  be  further  established  by  considering  the  out¬ 
age  range  II  criterion  in  light  of  the  measured  results  of  Fig.  5.1 
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and  th*  theoretical  results  of  TR12-76.  The  troposcatter  link  out 
age  duration  and  the  distribution  of  link  outage  duration  derived 
in  TR12-76  are  presented  here  as  Figs.  5.2.  and  5.3,  respectively. 
If  the  outage  range  II  criterion  is  to  be  met,  the  averaging  pro¬ 
cess  to  find  the  yearly  outage  probability  results  in  values  close 
to  hourly  outage  probabilities  for  a  Yq  °f  10dB.  At  lOdB  when  the 
r ms  Doppler  spread  (N)  is  5  Hz,  the  mean  outage  duration  is  less 
than  0.1  seconds.  Since  the  outage  duration  must  exceed  0.2 
seconds  to  count  in  range  II,  the  results  of  Fig.  5.3  show 
that  less  than  10%  of  the  outages  count  toward  an  outage. 

Thus  it  is  not  surprising  that  outage  range  II  was  not  found 
most  critical  when  5  Hz  was  used  as  the  C-band  rms  Doppler 
spread.  However,  the  median  value  of  measured  rms  Doppler 
spread  for  the  C-band  systems  of  Fig.  5.1  was  only  0.5Hz. 

A  crude  interpolation  of  the  mean  outage  duration  in  Fig.  5.2 
indicates  that  the  mean  would  be  significantly  greater  than 
0.2  seconds  and  hence  from  Fig.  5.3  we  can  expect  almost  all 
outages  tg.  }?e  in  outage  range  u. 

In  addition  to  these  factors,  the  limited  amount  of  rms 
Doppler  spread  data  makes  computation  of  outage  range  cri¬ 
teria  based  on  duration  statistics  subject  to  significant 
error.  For  these  reasons,  we  have  elected  a  conservative 
approach  which  uses  the  outage  range  II  criterion  for 
all  DCS  troposcatter  frequency  ranges  and  assumes 
that  all  outages  exceed  0.2  seconds  in  duration.  This  se¬ 
lection  is  only  slightly  conservative . for  C-band  systems 
because  few  outages  are  less  than  0.2  seconds  according  to 
the  data  presented  here.  Additionally,  frequent  recurring 
outages  of  the  type  in  ranges  IV  and  V  are  not  consistent 
with  the  slow  fading  measured  in  the  C-band  system  at  RADC. 


Fig.  5.2  Tropoacatter  Link  Outage  Duration 
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Using  this  criterion,  the  procedure  to  find  the  system 

gain  surplus  (deficit)  is  as  follows.  The  yearly  outage 

.4 

probability  for  a  BER  of  10  is  computed- for  diversity  con¬ 
figurations  2S/2F,  2S/2A,  and  2S/2F/2A.  In  the  long  term 
averaging,  this  calculation  uses  the  mean  multipath  spread 
for  the  quadruple  diversity  systems  and  an  assumption  of 

zero  multipath  spread  for  an  8th  order  diversity  system. 

"  •  •  ■  •*  •  ••  •  ••  m  -*•«  *  •  •  •  ,  . 

We  present  evidence  in  the  next  subsection  which  shows  that  the 
averaging  process  will  be  insensitive  to  the  2a  multipath  spread 
and  hence  little  error  is  introduced  by  using  the  mean  value 
for  the  quad  diversity  systems.  The  implicit  diversity  improve¬ 
ment  is  small  for  the  eighth  order  diversity  system,  and  hence 
the  use  of  zero  multipath  spread  is  slightly  conservative.  When 
the  yearly  outage  probability  is  plotted  on  the  short  term  BER 

distribution  as  a  function  of  y  ,  the  dB  shift  relative  to  the 

-4  ° 

yearly  objective  of  7.5  x  10  is  the  approximate  system  gain 
surplus  or  deficit.  The  result  of  this  procedure  is  presented 
in  the  next  subsection  for  the  RADC  test  link. 

b.  Computer  Program  Results 

The  method  of  computing  the  bit  error  rate  (BER)  statistics 
and  path  loss  distributions  for  a  troposcatter  equalizer  has 
been  applied  to  a  troposcatter  system  example.  The  path  para¬ 
meters  for  this  system  are  given  in  Table  5.2.  The  tropo- 
acatter  equalizer  chosen  for  the  analysis  is  the  MD-918  modem 
developed  by  Sylvania/SIGNATRON  for  U.S.  Army  ECOm  under  the 
Megabit  Digital  Troposcatter  Subsystem  [5.1].  The  modem  uses 
differentially  encoded  4PSK  modulation  with  a  four  channel 
adaptive  decision- feedback  equalizer  [5.2]  for  reception. 
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TH  P. 


Path  Distance  (mi.) 


Transmitter  Haight (ft.) 


Receiver  Haight (ft.) 


Transmitter  Antenna (ft.) 


Receiver  Antenna (ft.) 


Transmitter 

Boresight/Horizontal  <  (°) 
Boresight/Horizon  $  (°) 


Receiver 

Boresight/Horizontal  4  (°) 
Boresight/Horizon  4  (°) 


Squint  Angle  (°) 


Frequency  (6Hz) 


RADC 

Test-  Link 


168.3 


340. 


460. 


1.27 


0.27 


0.77 


0.27 


0.54 


Transmitter  power  -  kw 
N.F. 


TABLE  5 .3 


2a  MULTIPATH  SPREAD 


Squint 

Loss 


Main 


Elev. 

Beam 


Beamwidth 


RADC  (4.4-5.0GHZ) 


174. 


There  are  three  taps  per  forward  filter  diversity  channel  with 
tap  spacing  equal  to  one-half  a  4PSK  symbol  interval.  The 
backward  filter  uses  thre  ’  taps  in  its  transversal  filter 
realization  with  tap  spacing  equal  to  the  symbol  interval. 

Before  presenting  the  BER  statistics  for  our  example  link, 
the  short  term  BER  characteristics  of  the  DFE  modem  are  sunmarized. 
The  predicted  average  BER  performance  of  the  MO- 918  modem  as  a 
function  of  multipath  dispersion  is  shown  in  Fig.  5.4.  The 

■  •  •  i  . 

quantity  2<j/T  represents  the  ratio  of  twice  the  rms  multipath 
spread  to  the  4PSK  symbol  interval.  Mote  how  performance 
improves  with  increasing  2o/T  due  to  implicit  diversity  utili¬ 
zation  by  the  equalizer  but  then  performance  falls  off  again 
as  the  intersymbol  interference  penalty  exceeds  the  implicit 
diversity  gain.  The  amount  of  implicit  diversity  improvement 
and  reduction  of  intersymbol  interference  penalty  is  a  func¬ 
tion  of  the  number  and  spacing  of  forward  filter  taps  in  the 
DFE.  The  MD-918  uses  3  taps  with  a  spacing  of  T/2. 

The  BER  distribution  can  be  derived  in  closed  form  [5.3] 
when  the  rms  multipath  spread  is  zero,  i.e.  when  the  fading 

channel  is  non  frequency  selective  (flat  fading).  The  results 
of  this  computation  for  4PSK  transmission  and  a  Gaussian  noise 
detection  characteristic  comparable  to  DPSK  detection  are  shown 
in  Figs.  5.5  and  5.6  for  the  flat  fading  channel.  Ej^/Ng  is  a 
signal- to- noise  ratio  measure  defined  as  the  hourly  average 
received  energy  per  bit  divided  by  the  noise  power  in  a  1  Hz 
RF  bandwidth.  In  Section  4,  we  derived  for  the  first 
time  an  approach  for  the  calculation  of  the  short-term  BER 
distribution  when  the  channel  is  frequency  selective  fading. 
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Pig.  5.4  Short-T^rm  Av«rag*  BER 
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FLAT  FADING 
QUAD  DIVERSITY 


K)’5  10'6  10*7 


10'8  10'9 


ROR  RATE  (BER) 


Short -t -no  BER  Distribution 
Flat  Fading 


PROBABILITY  THAT  BER  IS  EXCE 


The  results  from  that  analysis  for  independent  equal  strength 
diversity  channels  are  shown  in  Fig*  5.7.  Note  that  in  the 
critical  probability  range  of  10  *  to  10  *  the  DFE  modem 
characteristic  is  not  very  sensitive  to  variations  in  the  rms 
multipath  spread.  This  result  allows  us  to  estimate  the  yearly 
outage  probability  accurately  without  knowing  the  form  of  the 
long  term  multipath  distribution. 

Before  tne  bit  error  rate  statistics  for  the  example  link 
can  be  computed*  squint  loss  is  calculated  and  the  2<j  multi¬ 
path  spread  is  determined  for  the  main  beam  and  elevated  beam* 
multipath  profiles.  These  values  are  given  in  Table  5.3  . 

The  data  rate  chosen  for  the  calculation  was  6.3  Mb/s  as 
it  corresponds  to  typical  user  requirements  for  strategic 
digital  troposcatter.  With  a  4PSK  modulation  format,  the  bit 
symbol  interval  T  is  then  317  nanoseconds.  The  2oq/T  value 
for  the  main  beam  diversity  is  0.42  for  the  RADC  path. 


The  results  of  the  bit  error  rate  computation  have  been  calculated 
as  a  function  of  the  mean  value  E^/N Q  for  a  2S/2F  diversity 
receiver.  To  accomplish  this  comparison  we  first  determine 
the  distribution  of  the  median  path  loss  using  the  common 
volume  integration  model  of  Section  3  and  the  NBS  time  avail¬ 
ability  analysis  discussed  in  Section  6.  The  common  volume 
integration  model  includes  the  effect  of  aperture  to  medium 
coupling  loss  and  provides  a  path  loss  parameter  normalized 
to  unit  gain  antennas.  To  convert  median  path  loss  into 
in  dB*  the  relationship  is 


i 


1 


4 
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Eb/No  PER  DIVERSITY,  dB 

Fig*  5*7  Short  T«rm  BER  Distribution 

for  Frequency  Selective  Fading 


VVVL-NF-10  log  (Rb)+i74+1«6 


VNo- 


where 

P_  ■  transmit  power  in  dBm 
T 

6T/R  "  antenna  gain  of  transmitter  receiver  in  dB 
L  ■  total  path  loss  including  coupling  loss  effects 
NF  “  noise  figure 

R^  “  data  rate  in  b/s 

174  *  noise  power  in  1  Hz  bandwidth  in  dBm 

1.6  ■  dB  amount  that  the  mean  exceeds  the  median  for 

a  Rayleigh  fading  channel. 


The  antenna  gain  formula  for  troposcatter  reflectors  of 
diameter  D  is 

G  -  10  log  (  tt2D2«/x2 

where  0  <  t  <  1  i>  the  antenna  efficiency  with  typical  value 
on  the  order  of  0.57. 

in  the  calculation  of  median  path  loss  the  refractive 
index  variance  dependence  is  based  on  Fried's  £5.3]  results 
at  optical  wavelengths.  Also  a  service  probability  of  0.95 
has  been  selected  as  a  conservative  measure.  We  have  not 

collected  sufficient  empirical  data  to  establish  the  absolute 

. . . 

accuracy  of  the  prediction  model.  For  the  RADC  test  link 
where  empirical  data  is  available*  the  winter  median  path 
losses  for  the  two  space  diversity  antenna  receivers  were 
measured  [5.5]  to  be  258  and  260  dB.  Our  angle  diversity 
test  program  from  October  to  May  1979  measured  an  average  of  the 
four  receivers  of  256  dB.  Our  predicted  value  for  this  link 
for  all  year  conditions  is  254  dB  which  is  good  agreement. 


1 


v 
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For  an  angle  diversity  application,  the  squint  loss  and 
the  availability  of  a  redundant  power  amplifier  has  been  con¬ 
sidered  in  the  system  comparison  with  conventional  systems. 

Since  the  2S/2F  configuration  requires  two  power  amplifiers, 
any  comparison  with  angle  diversity  should  be  on  the  same 
basis.  It  also  follows  that  a  2S/2A  system  would  require  a 
second  power  amplifier  for  failure  and  maintenance  redun¬ 
dancy  reasons.  For  these  practical  considerations  the  2S/2A 
system  utilizes  two  power  amplifiers  tuned  to  the  same  frequency.* 
Thus,  in  effect,  the  transmit  power  per  diversity  of  a  2S/2A 
system  is  3  dB  more  than  in  an  analogous  2S/2F  system.  The 
3  dB  advantage  of  the  2S/2A  system  over  the  2S/2F  system  is 
included  in  the  comparison  of  bit  error  rate  statistics. 

The  link  performance  is  summarized  in  a  graph  of  the  short 
term  BER  distribution  for  three  diversity  configurations  at  6.3  Mb/s: 
dual  space/dual  frequency  (2S/2F) ,  dual  space/dual  angle  (2S/2A) ,  and 
dual  space/dual  frequency /dual  angle  (2S/2F/2A) .  On  this  graph 
the  yearly  objective  for  the  outage  probability ,  the  yearly 
median  E^/Ng,  and  the  yearly  outage  probability  for  each  diver¬ 
sity  configuration  are  noted.  The  dB  shift  to  the  right  (left) 
of  the  yearly  outage  probability  relative  to  the  yearly  objective 
is  the  system  gain  surplus  (deficit).  Fig.  5.8  shows  these  results 
for  the  RADC  test  link.  Note  that  all  three  diversity  configura¬ 
tions  fail  to  meet  the  TR12-76  objective.  For  a  2S/2A  configura¬ 
tion  the  system  gain  deficit  is  approximately  3  dB.  Use  of  the 
eighth  order  system  reduces  this  deficit  to  about  2  dB.  Measured 

*Each  power  amplifier  drives  one  of the  two  transmit  antennas 
producing  a  composite  antenna  pattern  with  nulls  spaced  at 
intervals  equal  to  the  antenna  spacing.  Since  the  null  width 
is  small  compared  to  the  common  volume  dimension  defined  by 
the  beamwidth,  the  nulls  can  be  neglected.  Experimental 
tests  have  verified  the  usefulness  of  this  configuration. 


probability  that  be-r  is  exceeDE’o 


2S/2F  g^No  (<*B)  PER  DIVERSITY 
Fig.  5.8  RADC  Test  Link  Predicted  Performance 
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BER  statistics  summarized  in  the  next  section  illustrate  the 
negative  performance  margin  of  this  test  link  relative  to  the 
DCS  troposcatter  link  standards. 
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SECTION  6 

LONG  TERM  VARIABILITY 

Tropospheric  scatter  systems  are  subject  to  two  fading  phenom¬ 
ena  -  short  term  multipath  fading  and  long  term  power  fading.  The 
short  term  fading  of  the  instantaneous  received  power  within  periods 
of  time  ranging  from  less  than  a  second  to  many  minutes  results  from 
random  fluctuations  in  the  relative  phasing  between  component  waves 
arriving  at  the  receiver  over  slightly  different  propagation  paths. 
The  long  term  power  fading  results  from  slow  changes  in  average 
atmospheric  refraction,  in  the  intensity  of  refractive  index  turbu¬ 
lence,  and  in  the  degree  of  atmospheric  stratification.  The  power 
fading  is  characterized  by  hourly  or  diurnal  variations.  The 
evaluation  of  troposcatter  system  performance  is  accomplished  in 
part  by  determining  the  hourly  median  path  loss  where  the  median 
is  computed  to  include  the  short  term  multipath  fading  and  in  turn 
considering  the  median  path  loss  as  a  random  variable  subject 
to  a  power  fading  distribution.  The  median  path  loss  calculation 
utilizes  fixed  values  of  the  mean  and  variance  of  the  refractive 
index  and  an  assumption  on  the  degree  of  atmospheric  stratification. 
In  the  earlier  portion  of  this  report  we  have  presented  a  predic¬ 
tion  method  for  the  computation  of  median  path  loss.  This  method 
has  the  following  important  characteristics: 

•  The  average  atmospheric  refraction  is  fixed  by 
utilizing  an  effective  earth's  radius  of  K  times 
the  actual  radius  to  account  for  the  mean  refrac¬ 
tive  index. 

•  The  intensity  of  refractive  index  turbulence  is 
fixed  by  the  refractive  index  variance  which  is 
chosen  to  correspond  to  dry  winter  afternoons.* 

This  period  of  time  generally  experiences  the 
poorest  propagation  conditions. 


•This  time  period  is  referred  to  as  Time  Block  2  (TB2) 
in  the  NBS  prediction  method. 


•  The  atmospheric  structure  is  derived  from  turbu¬ 
lent  scattering  theory  which  leads  to  a  refractive 
index  spectrum  slope  of  m*ll/3.  This  structure 

is  more  applicable  for  higher  frequency  (>  1  GHz) 
troposcatter  systems  than  the  stratified  layer 
assumption  used  in  the  NBS  prediction  method  (m-5) . 

•  The  aperture-to-medium  coupling  loss  is  included 
as  an  integral  part  of  the  path  loss  calculations. 


For  this  method  the  power  fading  of  the  median  path  loss 

is  determined  by  variations  in  the  effective  earth  radius 

2 

factor  K,  variations  in  the  refractive  index  variance  a  ,  and 

n 

changes  in  the  atmospheric  structure  leading  to  other  values 

of  the  refractive  index  spectrum  slope  m.  Given  probability 

2 

density  functions  on  the  parameters  K,  »  and  m,  the  computa¬ 
tion  of  the  median  path  loss  long  term  distribution  would  be 
straightforward.  Unfortunately  there  is  little  empirical  data 
available  to  derive  such  densities.  Same  experimental  evidence 

and  analysis  TGossard*  1977]  indicate  that  the  refractive  index 
2 

variance  a  is  the  dominating  factor  in  producing  significant 
n 

variations  in  the  median  path  loss.  Development  of  experimental 
2 

data  on  c  over  long  periods  of  time  for  different  geographical 
n 

areas  would  provide  a  basis  for  predicting  long  term  variability 

of  troposcatter  systems.  At  the  present  time  the  only  method  of 

predicting  this  variability  is  to  use  path  loss  data  taken  from 

existing  systems  and  integrated  into  the  NBS  variability  model 

[Rice,  1967].  Much  of  the  empirical  path  loss  data  has  been  taken 

from  systems  with  operating  carrier  frequencies  below  1  GHz.  The 

performance  prediction  for  new  troposcatter  systems  operating 

in  the  4  to  5  GHz  frequency  region  may  be  subject  to  large 

errors  as  a  result.  However,  in  the  absence  of  empirical  data 
2 

on  either  a  or  median  path  loss  at  these  new  frequencies,  an 
n 

extrapolation  of  the  NBS  variability  model  is  the  only  realistic 


J 
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engineering  choice.  In  the  first  subsection  we  briefly  review 
the  NBS  model  and  describe  the  specific  parameters  integrated  into 
the  SIGNATRON  prediction  computer  program  for  computation  of 
long  term  variability.  In  the  subsequent  subsection  the  spe¬ 
cial  issues  which  must  be  addressed  in  the  evaluation  of  the 
long  term  variability  of  an  angle  diversity  system  are  considered. 

a.  NBS  Long  Term  Variability  Model 

Considerable  experimental  evidence  suggests  that  the  long 
term  distribution  of  the  hourly  median  path  loss  is  normally 
distributed  in  dB.  If  we  denote  L(p)  as  the  hourly  median  path 
loss  in  dB,  which  is  not  exceeded  p%  of  the  time,  the  normal 
distribution  is  defined  as  follows: 


median  path  loss 


L(p) 


J  exp  [(?“L<5°))2  /  2oL2  ] 


(6.1) 


where  L(50)  is  the  median  path  loss  computed  using  the  numerical 
integration  of  (3.2)  3.2  is  received  power  and  aL  is  the  long  term 
standard  deviation,  nps  Long  Term  Variability  Model  [Rice,  19671  uses  erpirical 

data  to  determine  L(10)  and  L(90)  from  which  normal  probability 
graph  paper  can  be  used  to  plot  L(p)  at  other  values  of  p. 

This  calculation  includes  the  effect  of  prediction  uncertainty 
through  a  parameter  called  service  probability.  The  service 
probability  is  the  probability  that  a  new  system  will  meet  the 
long  term  performance  predictions. 
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A  parameter  called  effective  distance,  dg,  has  been  found 

to  be  superior  to  other  parameters  such  as  path  length,  angular 

distance,  and  distance  between  actual  or  theoretical  horizons, 

in  predicting  the  long  term  variability.  The  effective 

distance  is  defined  in  terms  of  d_  which  is  the  distance  for 

81 

equal  diffraction  and  scatter  path  loss  and  dg  which  is  the 

smooth  earth  distance  between  radio  horizons.  These  parameters 
are  computed  from  the  following  relations 


d  -  65(100/f)1/3  tan 

al 


d  -  3  -3  </  2  hr 


(6.2) 

(6.3) 


where  the  effective  antenna  heights  lu  and  h_  are  expressed 

e  e 

in  meters,  the  path  length  d  is  in  kilometers,  and  the  radio 
frequency  f  is  in  MHz.  The  effective  distance  is  defined  in 
the  MBS  model  as 


130/  [1+  (d  -  d  )  /d]  km,  d 

*1  8o  so 

]  130  +  d_  -  d  km,  d_ 

^  -o  *1  *o 


<  d 


(6.4) 


A  variability  function  Y(p)  which  depends  on  the  effec¬ 
tive  distance  and  can  be  corrected  for  frequency  and  climate 
effects  is  used  to  determine  L(10)  and  L(90)  from  the  relation 


L(p)  -  L (50)  -  Y(p). 


(6.5) 


The  variability  function  is  defined  by 


Y  (p)  -  y0(p#  <ae)  g(p#  f) 

where  g(p,  f)  represents  an  average  of  many  effects  that  are 
frequency  sensitive.  The  MBS  variability  model  provides  empiri¬ 
cal  curves  and  an  analytic  function  representation  for  Yq (p,d#) . 
In  addition. the  model  uses  a  parameter  V(SO,d0)  to  adjust  the 
long  term  reference  median  path  loss  computed  by  the  NBS  method 
to  the  median  loss  L(50)  in  (6.5).  For  dry  winter  afternoons 
this  adjustment  is  0  dB. 

The  analytic  function  representation  for  YQ(p,de)  is  given 
for  dry  winter  afternoons  (Time  Block  2)  as 


V10-  V 

-vo<90.  a.) 


W]  «p(-c3d.  3)  +  VV 


(6.6) 


where 


W 


(1-f^/fJ  exp  (-c2 


(6.7) 


and  the  constants  have  values  (d  in  kilometers)* 

O 


♦The  values  for  f  and  fw  given  in  Annex  III,  Tech  Note  101  do 
not  agree  with  th®  continental  temperate  climate  variability  curve 
given  in  Section  10.5  of  Tech  Note  101.  The  quantities  f  and  fw 
in  the  table  have  been  adjusted  to  provide  good  agreement? 


Y  (10,  d  ) 
o  e 

-Y  (90,  d^ 

o  e 

C1 

1.04  x  10“5 

1.05  x  10“ 

C2 

4.28  x  10“8 

7.00  x  10“ 

C3 

3.51  x  10“8 

7.64  x  10“ 

nl 

2.71 

2.59 

a2 

2.91 

4.80 

n3 

3.41 

3.68 

f 

10.5 

8.5 

m 

f- 

3.0 

2.3 

For  frequencies  greater  than  400  MHz,  the  curves  for 
g(p,  f ) ,  winter  afternoons,  provided  in  the  MBS  model  can  be 
well  approximated  by 


g(10,  f)  *  g  (90,  f) 


f  1-0.6  LOG  (.0005  f)  400  <  f  <  2000 

\  1.0  2000  <  f . 

(6.8) 


For  a  normal  distribution,  other  points  on  the  loss  distri¬ 
bution  are  calculated  from 


Y(0.01)  * 

Y(0.1)  » 

Y(l.O)  » 


3.33Y(10) 

2.73YU0) 

2.00Y(10) 


Y(99.99) 

Y  (99.  9) 

Y  (99) 


2.90Y  (90) 
2.41Y  (90) 
1.82Y  (90)  • 


These  calculations  and  the  loss  distribution  resulting  from 
(6.5)  provide  the  long  term  variability  of  the -median  path  loss 
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for  a  continental  temperate  climate  during  winter  afternoons. 
The  UBS  variability  model  provides  numerous  curves  and  other 
variability  functions  to  accommodate  other  portions  of  the 
season  and  other  climatic  areas. 

The  path  loss  distribution  (6.5).  corresponds  to  a  service 
probability  of  50%,i.e.,  50%  of  the  systems  built  would  exceed 
the  performance  predictions.  Conservative  engineering  practice 
would  recommend  the  selection  of  a  service  probability  of  95%. 
For  this  choice  the  loss  distribution  must  be  adjusted  by  the 
prediction  error  according  to  the  NBS  formula 

L0  95<p)  -  Up)  +  1.65  \/ 12.73  +  0.12Y2 (p)  •  (6.9) 


The  variability  predictions  used  in  this  study  utilize  the 
NBS  long  term  variability  model  with  a  service  probability  of 
95%.  The  median  path  loss  L(50)  is  computed  from  the  previously 
described  numerical  integration  procedure  which  corresponds  to 
the  loss  for  turbulent  scattering  conditions  during  dry  winter 
afternoons  in  a  temperate  continental  climate.  The  formulas 
(6.1)  thru  (6.9)  are  then  used  to  derive  the  loss  distribution. 
The  mean  Eb/NQ  distribution  for  evaluation  long  term  variations 
in  digital  system  performance  is  determined  from  the  path  loss 
distribution  by  the  formula 


J*<P>  -  PT+  Gt+  Gr-  L>g5(p)  -[lOLOG  (1^)  +  NF  -  174]+ 1,6 


(6.10) 


where 


■  transmit  power  in  dBm 

■  transmit  antenna  gain  in  dB 
-  receiver  antenna  gain  in  dB 


ft 
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*  path  loss  not  exceeding  p%  of  the  time  for 
unit  gain  antennas  and  a  service  probability 
of  0.95 

■  data  rate  in  b/s 

■  receiver  noise  figure  in  dB 

■  received  noise  power  in  1  Hz  bandwidth 

*  factor  relating  median  to  mean  for  a  complex 
Gaussian  scatter  channel. 

0.  Long  Term  Variability  in  an  Angle  Diversity  System 

In  space  and  frequency  diversity  systems,  the  common  scat¬ 
tering  volume  is  virtually  the  same  for  each  diversity  and  hence 
the  long  term  median  path  loss  varies  the  same  for  each  diversity. 
In  angle  diversity  system,  however,  the  common  volumes  are  sepa¬ 
rated  and  the  long  term  variability  is  not  identical  for  the 
angle  diversity  beams.  The  effect  of  this  decorrelation  of  the 
diversity  power  fading  may  improve  the  system  availability  be¬ 
cause  a  power  fade  in  the  main  beam  diversity  is  not  always  ac¬ 
companied  by  a  power  fade  in  the  squinted  diversity  beam.Phys- 
ically  one  can  imagine  this  situation  in  a  vertical  angle 
diversity  system  where  the  inhomogeneous  structure  of  the  at¬ 
mosphere  results  in  say  a  larger  refractive  index  variance  or 
superior  atmospheric  stratification  at  heights  corresponding 
to  the  elevated  beam  than  at  heights  defined  by  the  main-beam 
common  volume.  Since  the  common  volumes  in  a  vertical  angle 
diversity  system  are  separated  by  approximately  one  beamwidth 
Q  at  a  distance  d/2  from  the  link  terminals, this  decorrelation 
results  from  atmospheric  variations  over  distances  on  the  order 
of  fld/2.  For  typical  troposcatter  applications  this  distance 
is  on  the  order  of  one  or  more  miles  at  a  height  above  the 


*  Aperture  to  medium  coupling  loss  of  the  actual  diameter 
antennas  is  included  in  this  parameter. 


l.95(p) 

*b 

NF 

-174dBm 
1.6  dB 
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earth's  surface  also  on  the  order  of  a  few  miles.  Experimental 
evidence  verifying  the  long  term  power  fading  decorrelation  in 
angle  diversity  systems  has  been  reported  on  by  [Monsen,  1972]  and 
[Troitskiy,  1972].'  Data  from  these  angle  diversity  experiments 
and  data  from  the  AAC  field  tests  will  be  used  to  establish  the 
availability  improvement  from  power  fading  decorrelation. 

The  long  term  medial  path  loss  or  received  signal  level  (RSL) 
can  be  described  by  a  multi-variated  normal  distribution  when  each 
path  loss  or  RSL  is  given  in  dB.  Path  loss  and  RSL  are  related  by 
a  sign  change  and  an  additive  constant.  We  use  RSL  in  the  subse¬ 
quent  analysis  because  emperical  data  is  more  often  expressed  in 
this  quantity.  We  consider  a  dual  vertical  angle  diversity  system 
with  or  without  additional  frequency  or  space  diversity.  Let 
represent  the  medial  RSL  in  dB  for  diversity  I  and  we  number  the 
diversities  such  that  the  main  beam  diversities  are  odd  and  the 
elevated  beam  diversities  are  even,  e.g.,  for  a  2S/2F/2A  system  we 
have 

x^,  x3  mainbeam  space  diversity  channels 

Xj,  x7  mainbeam  frequency  diversity  channels 

Xj,  x4  elevated  beam  space  diversity  channels  ~  . 

xc,  x0  elevated  beam  frequency  diversity  channels, 
o  o 

Because  of  the  power  fading  correlation  in  space  and  fre¬ 
quency  diversity,  one  has 

*1  "  X3  “  X5  *  *7 


X2  “  *4  *  X6  “  X8 


and  x^  and  x^  are  described  by  the  joint  normal  density  function 


f(Xj#  x2) 


1  -x’ Rx  /  2c' 

2nclR|ls 


(6.11) 


where 


ra.  ■  mean  value  of  median  RSL 

for  diversity  1  (main  beam) 

m_  ■  mean  value  of  median  SSL  • 

for  diversity  2  (elevated  beam) 

9  “  ma inbeam  RSL  standard 

deviation 

ao  ■  elevated  beam  RSL  standard  9 

deviation 

’  p  *  power  fading  correlation  coefficient. 

The  NBS  variability  model  discussed  in  the  last  subsection  '** 
provides  a  method  for  determining  the  path  loss  (RSL)  standard  devia¬ 
tion  a*  Since  p*o  and  a“l  for  space  and  frequency  diversity 
systems*  determination  of  9  completely  specifies  the  long  term 
variability  for  these  systems.  In  angle  diversity  systems  the 
parameters  p  and  a  will  also  influence  the  long  term  variability 
of  the  effective  received  signal  power. 


In  order  to  assess  the  effect  of  decorrelated  power  &ding 
and  an  increase  in  the  elevated  beam  standard  deviation*  we  use 
the  mean  diversity  RSL  in  dB  as  a  system  measure  of  per¬ 
formance.  This  measure  corresponds  to  the  geometric  mean  of 
the  diversity  a ignal-to-noise  ratios  which  asymptotically  for 
large  signal- to-noise  ratio  governs  the  short  term  performance 
of  both  analog  and  digital  systems.  Thus  we  define  the  mean 
diversity  RSL 


x 


+  x2) 


(6.12) 


and  a  system  availability  function  as  x(p)«  where  x(p)  is  the 
mean  diversity  RSL  which  is  exceeded  100%  of  the  time,  viz.. 
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P  *  prob  {mean  diversity  RSL  £  x(p)>. 

Since  the  sum  o£  two  normal  random  variables  is  also  normal , 
we  have 


P  =«  —1 —  /  exp  1 5  -  x)2  /  2a2]dC  (6.13) 

<Jx/2rr  x  (p) 

where  the  long  term  mean  and  variance  are  equal  to 
*  ■  \  t*i  +  *2) 

ox2  *  (x  -  x)2  =  j  {(Xj^  -  m^)  +  (x2  -  m2)]2 

2 

=  (1  +  2pa  +  a^)  . 

For  convenience  we  define  the  function 

«a,e  ■  t  1  *  2f  *  ^  J1/2  (6.14) 


which  reduces  to  unity  for  frequency  diversity  systems,  a=l,  p=l. 
The  availability  for  the  angle  diversity  system  is 


1  -  Q 


x(p) 

-  (m1+m2)/2 

7° 

Vl+2pa+a2 

=  1  -  Q 


ml+m2> 
x(p)  -  -ij-  2 


~oT 


a,  p 


(6.15) 


where  Q(.)  is  the  normal  distribution  function 

1  «  _-v2/2 


Q(u) 


u 

/  e 


dv 


Equation  (6.15)  for  angle  diversity  should  be  contrasted  with 
the  result  for  frequency  diversity  when  p«l,  a*l. 
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and  but  both  are  reduced  by  3  dB,  The  availability  for 

frequency  diversity  is 

x(p)  -  m.  +3 

P  =  1  -  Q  - ~= -  (6.16) 

The  reduction  of  the  mean  RSL  in  a  frequency  diversity  system 
relative  to  an  angle  diversity  system  results  from  the  necessity 
to  return  one  of  the  two  power  amplifiers  in  an  angle  diversity 
system  to  adjacent  frequency  for  a  frequency  diversity  application. 
This  action  reduces  the  received  power  for  diversity  by  3  dB. 

Let's  compare  the  frequency  (2S/2F)  and  angle  (2S/2A)  diversity 
systems  as  a  function  of  availability  level.  An  availability  of  p  is 
realized  in  the  two  systems  when  the  argument  of  Q  in  (6.15)  and 
(6.16)  is  equal  to  -qp  where 

P  -  1  -  Q(-qp) •  (6.17) 

For  equal  availabilities,  the  availability  levels  for  RSL  for  the 
two  systems  are  then 

m.  +m- 

XA(p)  =  — ^ —  "  qpafa,p  •  Angle  Diversity 

Xp(p)  *  m1  -  3-qpa  ;  Frequency  Diversity 

The  availability  improvement  in  dB  for  the  angle  diversity 
system  is 


The  first  term  is  the  power  amplifier  differential,  the  second 
is  the  squint  loss,  and  the  third  term  represents  the  combined  ef¬ 
fects  of  increased  standard  deviation  in  the  elevated  beam  (a>l) 

and  long  term  decorrelation  of  the  main  and  elevated  beam  common 
volumes  (p  <  1)  . 

In  (6.14)  one  can  see  that  the  effect  of  the  power  fading 

correlation  coefficient  is  to  improve  system  availability  with 

decreasing  values  of  the  coefficient  but,  on  the  other  hand,  an 
increase  (a  >  1)  in  the  elevated  beam  path  loss  standard  devia¬ 
tion  decreases  the  system  availability.  In  order  to  evaluate 
the  result  of  these  competing  effects,  we  examine  some  previous 
angle  diversity  experimental  data. 

Long  term  variability  data  was  taken  on  the  Bell  Laboratories 
experimental  angle  diversity  link  reported  on  by  [Monsen,  1972] 

The  important  characteristics  of  this  link  were 

•  triple  vertical  angle  diversity 

•  2.17  GHz  operating  frequency 

•  179  statute  mile  path 

•  28*  transmit  antenna,  50'  steerable  receive  antenna 

•  1  Kw  transmit  power 

•  1.6  beamwidth  separation  between  beams. 

Measurements  of  hourly  median  signals  once  or  twice  a  week 
over  a  period  from  September  1970  through  February  1971  provided 


a  set  of  43  independent  hourly  median  samples.  The  long  term 
received  power  means,  standard  deviations,  and  correlation  co¬ 
efficient  for  the  lowest  two  vertical  beams  have  been  computed 
and  are  given  in  the  following  table  as  a  function  of  antenna 
pointing  angle  relative  to  the  horizon.  The  last  addition  in  the 
table  is  the  computed  angle  diversity  improvement  at  an  avail¬ 
ability  of  99.9%.  (qp  =  3.1). 

For  large  take-off  angles  when  the  squint  loss  (X^  -  X2) 
is  large,  the  angle  diversity  system  is  seen  to  be  inferior  to  a 
frequency  diversity  system.  As  the  take-off  angle  is  reduced,  a 
point  is  reached  wehre  the  angle  diversity  is  superior.  In  this 
experimental  system  at  a  take-off  angle  of  0.5°,  the  angle  diversity 
system  is  1/2  dB  poorer  at  99.9%  availability.  By  reducing  the 
take-off  angle  to  0.16°,  thus  reducing  the  mean  received  power  by 

(89.3  -  86.7)  +  (91.3  -  93.0)  *  0.9  dB 

one  realizes  a  3.2  dB.  improvement  in  the  angle  diversity  system 
for  a  net  gain  of  2.3  dB.  Thus  an  angle  diversity  system  with  a 
take-off  angle  of  0.16  is  2.3  dB  better  than  a  frequency  diversity 
system  with  a  take-off  angle  of  0.50.  In  this  experiment  system, 
no  attempt  was  made  to  minimize  the  squint  angle.  A  feedhorn  with 
a  squint  angle  of  approximately  1.6  beamwidths  which  produced  a 
squint  loss  of  6.3  dB  at  the  take-off  angle  of  0.5.  Using  the 
feedhorn  design  developed  under  this  program,  a  reduction  in  the 
squint  loss  to  2  or  3  dB  is  judged  reasonable.  In  an  optimized 
design  the  angle  diversity  improvement  relative  to  frequency 
diversity  is  judged  on  the  basis  of  this  data  to  be  about  3  dB. 


main  beam  excessively  blocked 


When  the  squint  loss  is  equal  to  6  dB  there  is  exact  cancel¬ 
lation  of  the  3  dB  power  amplifier  gain  with  the  mean  diversity 
loss  associated  with  a  squinted  dual  angle  diversity  system.  For 
operational  systems  the  net  gain  due  to  these  two  factors  will 
normally  be  positive  since  the  antenna  beams  can  be  placed  closer 
together  than  the  1.6  beamwidth  separation  in  the  BTL  experimental 
system.  The  new  gain  due  to  the  power  fading  correlation  coef¬ 
ficient  and  elevated  beam  standard  deviation  is  more  difficult  to 
assess  because  of  the  rather  limited  empirical  data.  [Troitskiy, 
1972]  has  presented  data  on  three  links  which  show  a  smaller  power 
fading  correlation  coefficient  than  found  from  the  BTL  data.  Table 
6.2  summarizes  Troitskiy' s  experimental  results.  The  power  fading 
correlation  coefficient  was  measured  between  envelope  powers  in 
watts  rather  than  dBm  and  is  designated  p.  In  the  2nd  Interim 
report  [Monsen,  Pari,  1978]  we  determined  that  p  is  overbounded 
by  p  but  the  bound  is  tight  and  to  a  good  approximation  p  =  p 

It  is  surprising  that  the  results  in  Table  6.2  show  a  smaller 
.power  fading  correlation  coefficient  when  these  links  have  the 
beams  closer  together  in  the  sky  than  the  beams  were  in  the 
BTL  experimental  link.  Also  it  is  unfortunate  that  Troitskiy 
did  not  report  on  data  from  which  the  elevated  beam  standard 
deviation  can  be  calculated. 


Table  6.2 

Experimental  Results,  Ref.  [Troitskiy,  1972] 


Path 

Length 

Miles 

Freq. 

GHz 

width 

Angle  Between 

Beams 

mr 

Scatter 

Angle 

mr 

Hours 

of 

Data 

P 

188 

B 

m 

11.6 

35 

500 

0.50 

174 

8.7 

33 

240 

0.56 

267 

D 

B 

6.9 

50 

, 

100- 

150 

0.679 
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However,  the  smaller  value  of  long  term  correlation  suggests 
larger  values  of  angle  diversity  improvement  than  calculated  from 
the  BTL  data.  .« 

The  final  set  of  experimental  data  examined  as  part  of  an 
angle  diversity  availability  evaluation  is  the  field  teat  data 
taken  during  the  AAC  program.  The  RADC  link  for  the  AAC  testa 
had  the  following  characteristics: 

•  vertical  angle  diversity  in  a  2S/2F/2A  configuration 

•  4.5  and  4.69  GHz  receive  frequencies 

•  168  statute  mile  path 

e  28*  antennas 

e  2  to  4  Kw  transmit  power 

•  1.3  beamwidth  separation  between  beams. 

The  measurement  of  medial  received  powers  covered  a  period 
from  September  1977  to  October  1979.  The  statistical  results 
from  that  test  period  are  given  in  Table  6.3. 

Note  that  the  RADC  link  also  has  less  correlation  between  beams  that  measured 
in  the  Bell  Laboratories  experiment  and  the  beams  are  closer  together.  Outside 
of  statistical  and  measurement  error  effects,  this  result  suggests  a  more  homogeneous 
tropospheric  structure  in  the  S-band  BTL  tests  than  in  the  C-band  RADC  tests.  Also 
the  elevated  beam  RSL  has  less  variation  than  the  main  beam  in  the  RADC  experiments 
whereas  it  had  more  variation  in  the  BTL  experiments. 

To  compute  the  angle  diversity  availability  improvement  from 
(6.18)  we  reduce  the  data  in  Table  6.3  to 

Mean  Mainbeam  RSL  =  -94.1  dBm  *  m^ 

Mean  Elevated  RSL  =  -97.4  dBm  «  m2 

Mainbeam  Std.  Dev.  *  7.5  dB  ■  o 

Elevated  Beam  St.  Dev.  ■  6.3  dB  *  acr ,  a  =  0.84 

Avg.  Correlation  Coeff.  -  0.71  *  p 

The  availability  improvement  for  angle  diversity  at  a  99.9% 
availability  is 


Table  6.3 

AAC  Field  Test  Data 


X1 

0 

x2"*l 

ao 

RCVR 

MEAN 

STD. 

■MEAN 

STD. 

RSL 

SQUINT 

SQUINT 

A/S/F 

(RADC) 

dBm 

dB 

dB 

dB 

000 

-93.4 

8.5 

3.3 

5.5 

001 

-96.7 

7.2 

1.7 

5.4 

010 

i 

vo 

o 

• 

00 

7.6 

5.8 

5.0 

Oil 

-95.6 

6.5 

2.5 

5.7 

100 

-96.8 

6.4 

- 

- 

101 

-98.3 

6.5 

- 

- 

no 

-96.6 

5.9 

- 

- 

111 

-98.1 

6.4 

6-18 

XAF  "  3  "  +  (3,1)  (7*5)  <°-1487> 

■  4.8  dB. 


The  improvement  of  the  angle  diversity  system  over  frequency 
diversity  is  also  supported  by  a  direct  computation  of  the  avail¬ 
ability  from  the  actual  RSL  data  taken  in  the  RADC  experiment,  in 
the  RADC  experiment,  data  was  taken  on  pairs  of  frequencies,  an¬ 
tennas,  and  feedhoms.  From  this  data  base,  three  hypothetical 
systems  can  be  examined.  These  are 

1.  Conventional  2S/2F  using  4.50  and  4.69  GHz 

2.  Angle  diversity:  2S/2A  using  4.50  GHz 

3.  Angle  diversity:  2S/2A  using  4.69  GHz 


The  composite  RSL  for  each  of  these  systems  is  an  accurate 
measure  of  the  short  term  performance.  For  the  8  receivers  de¬ 
fined  in  Table  6.3  (A  *  angle,  S  *  space,  F  «*  frequency),  the 
composite  RSL  is  defined  as 


(x1+x2+x3+x4) 

2S/2F 

<X0+X2+X4+V 

+  3 

2S/2A 

4.50 

GHz 

(x1+x3+x5+x7) 

+  3 

2S/2A 

4.69 

GHz 

where  3dB  has  been  added  to  the  angle  diversity  systems  to  accu¬ 
rately  reflect  the  two  power  amplifier  configuration  of  a  traffic 
bearing  system.  The  long  term  variation  of  the  composite  RSL  (x) 
for  the  three  hypothetical  systems  is  shown  in  Fig.  6.1.  At  low 
RSL  values  both  angle  diversity  systems  show  a  significant 
improvement.  For  example  at  95%  availability,  the  angle  diver¬ 
sity  systems  improve  on  the  conventional  system  by  about  1V2. 
and  4dB.  Note  that  the  weaker  angle  diversity  system  has  a  poorer 
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median  composite  RSL,  but,  due  to  the  long  term  decorrelation 
effect,  the  angle  diversity  system  has  a  larger  RSL  at  high 
availability  values. 

The  empirical  results  presented  in  this  subsection  clearly 
indicate  the  system  gain  superiority  of  angle  diversity  systems 
over  frequency  diversity  systems.  This  system  gain  advantage 
is  complemented  by  the  2  to  1  bandwidth  advantage  of  the  angle 
diversity  systems. 
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